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PREFACE 


Laboratory  investigations  of  the  properties  of  clay  shales  by  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  were  requested 
and  authorized  by  the  Office,  Chief  of  Engineers,  U.  S.  Army,  in  1965 
under  Engineering  Study  (ES)  529  and  later  in  FY  1968  under  ES  542. 

The  study  reported  herein  was  conducted  under  CWIS  31151  and  CWIS  31244, 
"Strength  and  Deformation  Properties  of  Clay  Shales,"  during  the  period 
from  1975  to  1981. 

The  initial  phase  of  the  study  consisted  of  both  laboratory  test¬ 
ing  by  the  WES  and  the  review  of  laboratory  test  results  from  other  labo¬ 
ratories  of  the  Corps  of  Engineers.  Report  1  of  this  series,  "Develop¬ 
ment  of  Classification  Indexes  for  Clay  Shales,"  which  was  published  in 
June  1971,  summarized  existing  testing  procedures  for  evaluating  the 
index  and  physical  properties  of  clay  shales  and  provided  a  basis  for 
the  adoption  of  standard  pretreatment  procedures  (i.e.,  undried,  air- 
dried,  and  blenderized)  for  grain-size  determinations  and  Atterberg 
limits  tests  on  clay  shales. 

The  second  phase  involved  classification  indexes,  mineralogy,  and 
residual  shear  strength  testing  of  a  number  of  different  clay  shales. 
Report  2  of  this  series,  "Residual  Shear  Strength  and  Classification 
Indexes  of  Clay  Shales,"  which  was  published  in  August  1974,  compares 
various  laboratory  procedures  and  equipment  used  for  determining  the 
residual  strength  of  clay  shales,  and  the  effects  of  procedures  devel¬ 
oped  in  Report  1  on  the  classification  indexes. 

The  third  phase  concerned  the  concepts  and  laboratory  evaluation 
of  (a)  temperature  change  as  it  affects  the  development  of  pore  pres¬ 
sures  in  clay  shale,  (b)  changes  in  pore  pressure  under  incremental 
isotropic  and  axial  stress  changes  in  triaxial  test  specimens,  and  (c)  ef 
fects  of  sample  anisotropy  on  pore  pressure  development.  Report  3  of 
this  series,  "Preliminary  Triaxial  Test  Program  on  Taylor  Shale  from 
Laneport  Dam,"  which  was  published  in  September  1976,  summarizes  the  tri¬ 
axial  testing  program  results,  the  problems  encountered  in  testing  clay 
shales,  and  the  concepts  used  in  analyzing  the  clay  shale  triaxial  test 


This  report,  phase  four  of  the  study,  consists  of  defining  the 
pore  pressure  and  volume  change  characteristics  of  clay  shales  and  devel 
oping  a  technique  for  computing  excess  construction  pore  pressures  in 
clay  shale  foundations  through  the  use  of  laboratory,  analytical,  and 
field  investigations.  The  laboratory  testing  was  conducted  by  Mr.  P.  A. 
Gilbert  of  the  Soils  Research  Facility  (SRF),  Soil  Mechanics  Division 
(SMD) ,  Geotechnical  Laboratory  (GL),  WES.  The  analysis  of  the  labora¬ 
tory  and  field  data  was  conducted  by  Messrs.  D.  A.  Leavell  and  J.  F. 
Peters,  SRF.  The  laboratory  testing  procedures  and  theoretical  concepts 
used  in  data  analysis  evolved  from  contributions  by  Dr.  R.  H.  G.  Parry 
in  the  third  phase  of  the  study.  The  computer  program  CURLS  and  the  in¬ 
fluence  chart  were  developed  by  Mr.  Peters.  The  field  data  for  Hills¬ 
dale  Dam  were  obtained  from  the  U.  S.  Army  Engineer  District,  Kansas 
City.  Mr.  Rollie  Fehrman,  Foundations  and  Materials  Branch,  Kansas  City 
District,  assisted  in  the  interpretation  of  these  field  data.  Mr.  Ralph 
R.  W.  Beene  was  technical  monitor  throughout  this  investigation.  From 
his  engineering  experience  in  the  behavior  of  clay  shales  he  has  con¬ 
tinually  provided  technical  guidance  for  this  study. 

This  report  was  prepared  by  Messrs.  Leavell  and  Peters  and 
Dr.  F.  C.  Townsend,  former  employee  of  the  WES,  under  the  general  super¬ 
vision  of  Mr.  C.  L.  McAnear,  Chief,  SMD;  Mr.  J.  P.  Sale,  former  Chief, 

GL  (retired);  and  Dr.  W.  F.  Marcuson  III,  Chief,  GL. 

COL  John  L.  Cannon,  CE,  COL  Nelson  P.  Conover,  CE,  and 
COL  Tilford  C.  Creel,  CE,  were  the  Commanders  and  Directors  of  the  WES 
during  the  conduct  of  this  investigation  and  preparation  of  this  report. 
Mr.  Fred  R.  Brown  was  the  Technical  Director. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


I 


’  i 

t 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 


verted  to  metric  (SI) 

units  as  follows: 

Multiply 

By 

To  Obtain 

degrees  (angular) 

0.01745329 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degree: 
Kelvins* 

feet 

0.3048 

metres 

inches 

25.4 

millimetres 

miles  (U.  S.  statute) 

1.609347 

kilometres 

pounds  (force)  per 
square  foot 

47.88026 

pascals 

pounds  (force)  per 
square  inch 

6894.757 

pascals 

pounds  (mass)  per 
cubic  foot 

16.01846 

kilograms  per 
cubic  metre 

<  • 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read¬ 
ings,  use  the  following  formula:  C  =  (5/9) (F  -  32).  To  obtain  Kelvin 
(K)  readings,  use:  K  =  (5/9) (F  -  32)  +  273.15. 

«  • 


•  • 
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ENGINEERING  PROPERTIES  OF  CLAY  SHALES 


LABORATORY  AND  COMPUTATIONAL  PROCEDURES 
FOR  PREDICTION  OF  PORE  PRESSURES 
IN  CLAY  SHALE  FOUNDATIONS 

PART  I :  INTRODUCTION 

Background 

1.  Design  and  stability  analyses  of  structures  founded  on  clay 
shales  require  knowledge  of  the  magnitude  of  pore  pressures  induced  by 
applied  loads  (Beene  1967,  and  Jasper  and  Peters  1979).  Many  clay 
shales  are  anisotropic  and  virtually  impervious.  These  two  factors  in 
combination  often  lead  to  the  development  of  high  pore  pressures  that 
do  not  readily  dissipate.  As  shown  in  Report  2  of  this  series  (Parry 
1976) ,  the  pore  pressure  response  of  clay  shales  can  be  predicted  from  the 
theory  of  transverse  isotropic  elasticity  by  assuming  the  clay  shale  to  be 
undrained.  Further,  preliminary  test  data  from  a  dam  of  moderate  height 
suggest  that  the  elastic  pore  pressure  response  is  valid  for  applied 
stresses  considerably  greater  than  those  imposed  by  the  embankment. 

2.  The  construction  failure  of  the  Waco  Dam  provides  an  excellent 
example  of  the  problems  associated  with  embankment  construction  on  clay 
shales.  The  failure  resulted  from  block  sliding  within  the  Pepper  shale 
formation,  which  caused  a  slide  within  the  embankment  (Beene  1967,  and 
Little  1968).  The  failure  zone  was  evidently  an  inherently  weak  horizon¬ 
tal  seam  within  the  shale  foundation.  From  piezometers  installed  after 
the  failure,  it  was  observed  that  excess  pore  water  pressures  were  on  the 
order  of  70  to  100  percent  of  the  added  embankment  load  (Beene  1967) .  Re¬ 
sults  of  stability  analyses  indicate  that  high  induced  pore  water  pressure 
contributed  to  this  slide  (Beene  1967,  and  Little  1968).  As  noted  by  Parry 
(1976) ,  the  pore  water  pressure  may  not  have  been  generated  within  the 
failure  zone  but  within  the  surrounding  intact  shale.  For  example, 

Beene  (1967)  reported  that  piezometers  located  outside  the  failure  area 


showed  the  pore  pressures  induced  beneath  the  stable  area  of  the  dam  were 
100  percent  of  the  embankment  load.  Therefore,  these  observed  pore  pres¬ 
sures  were  a  result  of  the  elastic  or  nearly  elastic  response  of  the 
clay  shale. 

3.  To  establish  a  procedure  to  predict  pore  pressure  response  in 
clay  shale  a  laboratory  investigation  was  performed  on  the  Taylor  shale 
from  Laneport  Dam,  Texas  (Parry  1976).  One  of  the  major  findings  in 
this  investigation  was  that  when  loaded  within  the  elastic  range,  this 
material  displayed  transverse  isotropic  elastic  behavior.  Further,  by 
relating  the  pore  pressure  response  to  the  elastic  behavior  of  the  clay 
shale,  the  predicted  pore  pressure  response  was  found  to  compare  favor¬ 
ably  with  that  measured  during  undrained  loading.  The  subsequent  work 
described  herein  has  been  directed  toward  developing  general  procedures 
for  evaluating  pore  pressure  response  in  clay  shales. 

Purpose  and  Scope 

4.  The  purpose  of  this  investigation  is  to  characterize  the  pore 
pressure  and  volume  change  response  of  clay  shales  and  develop  a  tech¬ 
nique  for  predicting  pore  pressure  response  in  the  field. 

5.  This  report  presents  theoretical  concepts  and  supporting  labo¬ 
ratory  and  field  data  necessary  for  relating  laboratory  behavior  to  mea¬ 
sured  field  pore  pressures  induced  by  changes  in  loading  conditions. 
Changes  in  loading  conditions  caused  by  excavation  and  embankment  place¬ 
ment  were  of  primary  interest;  however,  laboratory  loadings  up  to  fail¬ 
ure  are  also  presented. 

6.  To  verify  the  theoretical  concepts  and  obtain  background  data, 
four  different  clay  shales  were  subjected  to  cycles  of  loading  and  un¬ 
loading.  A  simplified  procedure  was  employed  to  compute  the  changes  in 
loading  caused  by  excavation  and  embankment  placement.  This  procedure 
combined  with  the  laboratory  pore  pressure  response  was  used  to  predict 
field  pore  pressure  behavior.  Theoretical  concepts  were  used  to  assess 
the  probable  accuracy  of  this  prediction  technique.  A  detailed  correla¬ 
tion  was  made  between  the  predicted  and  the  observed  pore  pressure 
response  using  piezometric  data  from  the  Hillsdale  Dam  project. 
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PART  II:  TESTING  PHILOSOPHY  AND  THEORETICAL  CONCEPTS 

General  Approach 

7.  The  method  used  to  evaluate  pore  pressure  response  is  based  on 
the  pore  pressure  parameter  concept.  A  pore  pressure  parameter  is  a 
numerical  constant  that  relates  the  change  in  pore  pressure  to  a  change 
in  applied  total  stress.  In  an  elastic  material,  the  complete  pore 
pressure  response  can  be  adequately  defined  by  two  parameters  (Skempton 
1954)  as  follows: 

Au  =  B[Ao3  +  ACAOj  -  Act3)]  (1) 

where 

Au  =  change  in  pore  pressure 
B  =  Skempton’ s  B  parameter 
A  =  Skempton' s  A  parameter 
ACTj,  Aa3  =  principal  stress  change  (total  stress) 

Both  the  A  and  B  parameters  can  be  measured  directly  in  a  triaxial 
compression  test,  but  it  is  important  to  note  that  under  the  assumption 
of  no  drainage,  they  can  also  be  related  to  the  elastic  properties  of 
the  material.  For  an  anisotropic  material,  such  as  clay  shale,  the  A 
parameter  depends  on  the  orientation  of  the  sample  with  respect  to  the 
applied  principal  stresses.  To  determine  the  A  parameter  in  the  labo¬ 
ratory  for  general  loading  conditions  numerous  laboratory  tests  would  be 
required.  However,  the  A  parameter  can  be  related  to  the  anisotropic 
elastic  constants.  An  appropriate  value  of  A  can  be  determined  for 
any  loading  case  encountered  in  the  field  provided  the  elastic  constants 
are  known.  Thus,  emphasis  has  been  given  to  laboratory  determination  of 
both  pore  pressure  parameters  and  anisotropic  elastic  constants. 

Elastic  Constant  Determination 


8.  The  philosophy  of  the  testing  program  was  to  subject 
men  to  a  variety  of  loadings  for  the  complete  determination  of 
tic  properties  under  triaxial  stress  conditions.  The  need  for 


a  speci- 
its  elas- 
data  from 


'  • 
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multiple-load  applications  can  best  be  understood  from  Figure  la,  which 
depicts  the  assumed  anisotropic  nature  of  the  clay  shales  tested.  In 
general,  five  independent  constants  are  needed  to  characterize  a  trans¬ 
versely  isotropic  material  completely  (Figure  la) .  For  the  typical 
field  problem  in  which  plane  strain  is  assumed  (Figure  lb),  only  four 
independent  constants  are  necessary  since  two  constants  can  be  expressed 
as  a  ratio.  However,  for  the  triaxial  text,  only  three  constants  can  be 
measured  independently  (Figure  lc),  since  two  of  the  applied  stresses 
are  equal.  As  shown  in  Appendix  A  (Equation  A13),  the  elastic  constants 
measured  in  the  triaxial  test  can  be  related  by 

Ae  =  Act  C  +  2 (Act  +  Act  )C  +  2Act  C  (2a) 

v  a  aa  a  r  ar  r  rr 


Ae  -  Ae  =  Act  C  +  (2Act  -  Act  )C  -  Act  C  (2b) 

a  r  a  aa  r  a  ar  r  rr 

where 

Ae  =  Ae  +  2Ae  =  change  in  volumetric  strain 
v  a  r  ° 

Aea  =  change  in  axial  strain 

Ae^  =  change  in  radial  strain 

Acr^  =  change  in  axial  effective  stress 

Act^  =  change  in  radial  effective  stress 

C  =  constant  relating  purely  axial  stress  to  axial 
aa  strain 


C  =  constant  relating  purely  axial  stress  to 
radial  strain 

C  =  constant  relating  purely  radial  stress  to 
radial  strain 


Equations  2a  and  2b  involve  three  constants;  thus  some  combinations  of 
three  equations  from  at  least  two  independent  loading  increments  are  re¬ 
quired  for  their  determination. 

9.  The  typical  loading  procedure  used  to  determine  the  elastic 
constants  is  depicted  by  the  stress  path  in  Figure  2.  The  sample  was 
first  subjected  to  isotropic  drained  loading  (initial  consolidation)  as 


DEVIATOR  STRESS 


TYPICAL  FIELD  PROBLEM 
(PLANE  STRAIN) 

LABORATORY  TEST 


FIVE  INDEPENDENT 
CONSTANTS 

(a) 


Figure  1.  Elastic  constants  required  to  characterize  clay  shale 
behavior  for  field  and  laboratory  loadings 


Figure  2.  Idealized  stress  paths  used  to  test 
clay  shale  specimens 


shown  by  increment  0.  The  specimen  was  next  loaded  axially  undrained  as 
indicated  by  increment  1.  Increment  2  represents  subsequent  consolida¬ 
tion  of  the  sample  as  shown  by  the  horizontal  effective  stress  path  and 
unchanging  total  stress.  An  undrained  increment  followed  by  an  incre¬ 
ment  of  consolidation  will  be  referred  to  as  a  loading  step  throughout 
this  report.  Numbers  3  through  8  along  the  path  represent  similar  in¬ 
crements  of  undrained  axial  loading  followed  by  consolidation  under 
constant  total  stress. 

10.  Each  of  the  undrained  and  drained  increments  of  the  stress 
path  yields  the  two  independent  equations  (2a  and  2b).  Thus,  each  set 
of  undrained  and  drained  increments  provided  sufficient  data  to  evaluate 
the  three  elastic  constants.  For  example,  using  Equation  2a  for  an  un¬ 
drained  increment  (Aev  =  0)  and  Equations  2a  and  2b  for  a  drained  incre¬ 
ment  results  in  three  independent  equations  with  three  unknown  constants. 
Following  this  procedure  for  each  pair  of  increments  in  the  stress  path 
gives  the  constants  for  each  step.  If  the  specimen  was  purely  elastic, 
the  constants  would  be  the  same  for  each  step.  However,  it  was  found 
that  the  elastic  constants  varied,  sometimes  systematically,  throughout 
the  cycle  of  loading  and  unloading.  Thus,  some  judgment  is  required  to 
determine  which  measured  value  is  most  representative  of  the  expected 
field  behavior. 

Pore  Pressure  Parameter  Determination 


11.  The  pore  pressure  parameters  were  obtained  from  the  pore 
pressure  induced  during  each  undrained  increment.  The  B  parameter  was 
generally  on  the  order  of  1.0;  thus  the  A  parameter  could  be  deter¬ 
mined  directly  from  Equation  1,  that  is 

Au  -  Act 

A  =  Ao  -  Aa  (3) 

a  r 


The  A  parameter  is  uniquely  defined  by  the  elastic  constants  and  should 
therefore  be  the  same  for  any  loading.  However,  as  in  the  case  of  the 
elastic  constants,  small  variations  in  the  A  parameter  were  observed 
for  each  step  of  the  loading  cycle. 
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PART  III:  DESCRIPTION  OF  TESTING  PROGRAM 

Equipment 

12.  Triaxial  tests  were  performed  in  the  stress-controlled  sys¬ 
tems  shown  in  operation  in  Figure  3.  The  specimen  was  loaded  axially 
from  the  top  by  a  double-acting  pneumatic  loading  system.  The  top 
platen  was  rigidly  connected  to  the  loading  rod  to  permit  loading  in 
either  compression  or  extension.  Cell  pressure  was  also  controlled 
pneumatically,  with  silicon  oil  used  as  a  confining  fluid  to  protect 
internal  instrumentation. 

13.  The  back  pressure  system  is  shown  schematically  in  Figure  4. 
Diffusion  of  compressed  air  into  the  back  pressure  system  in  tests  of 
long  duration  caused  degradation  of  the  B  value;  therefore,  an  air  bar¬ 
rier  was  necessary  to  isolate  the  back  pressure  system  and  the  saturation 
fluid.  The  air  barrier  for  the  saturation  reservoirs  consisted  of  an  im¬ 
pervious  steel  membrane  with  O-ring  seals.  The  burette  fluid  was  sepa¬ 
rated  from  the  air  by  a  1-  to  2-in.*  mercury  barrier.  The  mercury  bar¬ 
rier  was  separated  from  the  water  in  this  burette  by  dyed  kerosene  to 
facilitate  accurate  volume  change  measurements.  Also,  a  similar  air  bar¬ 
rier  was  used  to  isolate  the  confining  fluid  from  the  pneumatic  pressure 
system.  This  system  evolved  during  the  testing  program  and  all  these 
features  were  operational  only  during  the  testing  of  Pierre  and  Bearpaw 
shales . 

14.  Instrumentation  consisted  of  (a)  an  externally  mounted  load 
cell,  (b)  linear  variable  differential  transformer  (LVDT)  clamps,  lo¬ 
cated  at  specimen  one-third  points,  for  measuring  axial  and  radial  defor¬ 
mations,  (c)  an  external  LVDT  to  provide  axial  deformation  measurements, 
(d)  three  pore  pressure  transducers  (high  compliance)  monitoring  pore 
pressures  at  the  top,  bottom,  and  interior  of  the  specimen,  (e)  a  pres¬ 
sure  transducer  measuring  chamber  pressure,  and  (f)  a  digital  clock. 

Data  from  the  LVDT's  and  time  were  recorded  on  a  Digitrend  40-channel 

*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ments  to  metric  (SI)  units  is  presented  on  page  5. 
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Figure  3.  Laboratory  equipment  in  operation 


Figure  4.  Schematic  diagram  of  pore  water  pressure  (PWP) /drainage  system 


digital  printer.  Volume  changes  were  measured  during  consolidation  us¬ 
ing  a  small  bore  burette  read  by  a  vernier  cathodometer  to  the  nearest 
0.0026  cc. 

15.  Previous  work  on  clay  shales  (Parry  1976)  indicated  that 
changes  in  temperature  severely  influence  pore  pressure  measurements, 
therefore  testing  was  performed  in  an  environmental  room  (±1°F)  to 
minimize  temperature  effects. 

Materials  Tested 


16.  The  four  shales  used  in  this  study  and  their  source  locations 

are:  (a)  Bearpaw,  Billings,  Montana;  (b)  Kincaid,  Cooper  Dam,  North¬ 

eastern  Texas;  (c)  Pierre,  Limon,  Colorado;  and  (d)  Quivira,  Hillsdale 
Dam,  Kansas. 

17.  The  Bearpaw  and  Pierre  shales  were  sampled  for  the  Federal 
Highway  Administration  expansive  soil  investigation  described  by  Snethen 
(1979).  Their  source  locations  correspond  to  Sites  19  and  14,  respec¬ 
tively,  in  Snethen' s  report.  The  Kincaid  shale  was  obtained  from  a 
piezometer  boring  in  the  spillway  of  Cooper  Dam,  northeastern  Texas. 

The  Quivira  shale  was  sampled  from  the  cutoff  trench  of  Hillsdale  Dam, 
Kansas,  at  sta  99  +  50.  The  physical  characteristics  of  these  shales 
are  summarized  in  the  following  tabulation.  Each  shale  was  classified 
as  either  soft  or  stiff  based  on  a  visual  inspection.  Bearpaw  and 
Kincaid  shales  appeared  stiff,  whereas  Pierre  and  Quivira  shales  were 
comparatively  soft. 


Shale 

Bearpaw 

Kincaid 

Pierre 

Quivira 


Note:  The  definitions  of  the  characteristics  are:  LL  =  liquid  limit, 
PL  =  plastic  limit,  PI  -  plasticity  index,  Gs  =  specific  gravity 
of  solids,  wnat  =  natural  water  content,  Y<j  =  dry  unit  weight, 
and  eQ  *  initial  void  ratio. 

*  Average  values  (U.  S.  Army  Engineer  District,  New  Orleans  1977). 


w  Y 

Blenderized  _  nat  d 


LL 

PL 

PI 

U 

S 

percent 

pcf 

O 

58 

24 

34 

2.74 

16.6 

116.0 

wm 

85* 

65* 

2.72* 

19.5 

103.0 

mm 

58 

18 

40 

2.66 

22.9 

101.0 

mm 

51 

24 

27 

2.80 

14.8 

121.0 

0.44 

Test  Procedures 


18.  Preparation  of  triaxial  test  specimens  consisted  of  (a)  trim¬ 
ming  the  specimen,  (b)  inserting  the  interior  pore  pressure  probe, 

(c)  assembling  the  triaxial  chamber,  (d)  back  pressure  saturation, 

(e)  consolidation,  and  (f)  stress  path  loading. 

19.  The  specimens  were  removed  from  their  protective  shipping  con¬ 
tainers  and  roughly  trimmed  to  their  approximate  final  dimensions  using 

a  band  saw.  The  finished  2.8-in.  specimen  diameter  was  obtained  by  hand 
trimming  with  knives  and  straightedges  using  a  soil  trimming  lathe.  The 
specimen  end.;  were  trimmed  perpendicular  to  the  specimen  sides  in  a 
miter  box  to  the  final  6-in.  length.  All  final  trimming  was  performed 
in  a  humidity-controlled  room  to  preserve  the  specimen's  water  content. 

20.  The  pore  pressure  probe  was  fabricated  from  a  1/16-in.- 
outside-diam  stainless  steel  tube.  The  tip  of  the  tube  was  shaped  into 
a  point  and  covered  with  a  protective  coat  of  epoxy.  The  probe  tip  was 
then  finished  by  cutting  small  grooves  into  the  sides  of  the  tip.  To 
prevent  splitting,  the  specimen  was  confined  in  a  split  mold  before  the 
internal  pore  pressure  probe  was  inserted.  A  pilot  hole  was  drilled  in 
the  specimen  at  a  60-deg  angle,  to  within  1/2  in.  of  the  final  probe 
depth,  using  a  1/16-in.  standard  metal  drill  bit.  The  probe  was  care¬ 
fully  pressed  to  its  final  depth  of  1.5  in.  from  bottom  of  specimen. 

The  shank  of  the  probe  was  sealed  with  epoxy  to  the  base  platen  as  the 
specimen  was  mounted  in  the  cell. 

21.  The  specimens  were  set  up  by  a  procedure  similar  to  that  out¬ 
lined  for  R  tests  in  EM  1110-2-1906  (Dept,  of  the  Army,  Office,  Chief  of 
Engineers,  1970).  One  notable  addition  to  the  standard  procedure  was 

the  placement  of  LVDT  measurement  clamps  at  specimen  one-third  points. 

The  clamps  were  positioned  and  epoxied  to  the  membrane  after  a  vacuum 
had  been  applied  to  the  specimen.  The  standard  procedures  for  satura¬ 
tion  and  consolidation  were  used. 

Interpretation  of  Test  Data 

22.  The  parameters  measured  in  this  study  varied  with  time  as  the 


sample  reached  an  equilibrium  after  each  load  application.  Often  fluc¬ 
tuations  within  the  pressure  system  created  small  variations  in  the  ap¬ 
plied  loading  during  the  period  required  for  the  equilibrium  to  be 
obtained.  Therefore,  some  judgment  was  involved  in  selecting  the  best 
data  for  use  in  the  analyses. 

23.  The  test  data  from  a  load  increment  on  the  Pierre  Shale  give 
a  good  example  of  typical  fluctuations  in  measured  parameters  with  time 
(Figures  5,  6,  and  7).  Figure  5  shows  that  the  pore  water  pressure 
reached  an  equilibrium  within  2  min  at  the  top  and  bottom  of  the  speci¬ 
men  but  required  approximately  60  min  to  attain  an  equilibrium  at  the 
center  of  the  specimen.  Also,  even  with  fluctuations  in  the  axial  and 
radial  stresses,  the  pore  water  pressure  (Figure  5)  and  the  A  param¬ 
eter  (Figure  6)  were  not  significantly  affected.  The  strains  in  Fig¬ 
ure  7  did  not  reach  an  equilibrium  as  rapidly  as  the  pore  water  pressure 
and  A  parameter.  The  dashed  lines  in  Figures  5,  6,  and  7  denote  a 
loss  of  system  pressure.  Even  though  a  permanent  volumetric  strain  was 
induced  into  the  specimen  by  this  pressure  loss,  the  A  parameter  re¬ 
mained  in  the  0.7  to  0.8  range. 

24.  The  reliability  of  the  strain  data  was  also  of  considerable 
importance;  therefore,  several  redundant  measurements  were  used.  The 
axial  strain  was  measured  both  on  the  specimen  and  by  the  deflection  of 
the  loading  rod.  Further,  the  volumetric  strain  during  the  drained  load 
ing  step  could  be  determined  directly  from  burette  readings  or  computed 
from  the  radial  and  axial  strains.  Figure  8  presents  a  comparison  of 
the  volumetric  strain,  computed  from  axial  and  radial  strain,  and  the 
volumetric  strain  determined  from  the  burette  reading.  Note  that  the 
volumetric  strains  computed  from  linear  strain  measurements  agree  quite 
well  with  the  measured  volume  changes.  Further,  axial  strains  computed 
from  the  exterior  LVDT  data  agree  well  with  the  computed  strains  using 
the  interior  LVDT  data. 

25.  Where  inconsistencies  were  noted  in  redundant  measurements, 
the  following  order  of  preference  was  used: 


strain  (io“4) 
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Figure  7.  Strains  versus  time  after  stress  application 


Figure  8.  Burette  reading  versus  deformation  measurements  for 
volumetric  strain  calculations  during  drained  increments 
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Preferred  Order  of  Reliabilit' 


Axial  strain 


Volumetric  strain 


Radial  strain 


Pore  pressure 


Outside  LVDT 
Inside  LVDT  clamp 

Burette 

Computed  from  axial  and  radial 
strains 

Computed  front  burette  and  axial 
strain 

Inside  LVDT  clamp 

Bottom  stone 
Top  stone 
Probe 
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26.  Another  factor  requiring  some  judgment  was  the  selection  of 
the  loading  step  to  use  in  determining  the  elastic  constants.  Any  given 
increment  would  provide  only  two  of  the  three  independent  equations  re¬ 
quired  for  calculating  the  constants.  For  a  linear  elastic  material, 
there  is  no  such  difficulty  since  the  elastic  properties  would  be  the 
same  for  all  steps.  However,  for  the  soft  clay  shales,  the  repeated 
loadings  tended  to  induce  nonrecoverable  strains  (Figure  9).  Partic¬ 
ular  difficulty  was  encountered  when  attempting  to  incorporate  data 
from  the  initial  isotropic  consolidation  step  into  the  calculation 
process.  For  this  reason,  the  data  from  an  undrained  increment  were 
combined  with  the  data  from  the  subsequent  drained  increment  to  de¬ 
termine  all  three  constants  reliably.  Figure  10  shows  that  this  problem 
is  not  as  severe  in  the  stiff  shales. 


PART  IV:  TEST  RESULTS 


Determination  of  Elastic  Constants 


27.  The  data  for  the  triaxial  tests  on  the  four  clay  shales  are 
presented  in  Tables  1-4.  The  values  given  for  stress,  strain,  and  the 
A  parameter  represent  values  obtained  after  the  specimen  reached  an 
apparent  equilibrium.  The  axial  strain  data  and  pore  pressure  data, 
which  were  measured  at  various  locations  on  the  specimen  (see  para¬ 
graph  14),  are  the  values  considered  to  be  the  most  accurate  based  on 
the  order  of  preference  in  paragraph  25  and  are  not  averages  of  all 
values  measured. 

28.  The  procedure  used  to  compute  the  elastic  constants  can  best 
be  illustrated  by  performing  an  example  calculation.  In  Table  3,  the 
following  data  are  given  for  Pierre  shale: 

Increments  3  and  4 


Increment  3 
(Undrained) 


Equation  2a  for  the  undrained  quantities  and  Equations  2a  and  2b  for 
the  drained  quantities  give  the  following  three  independent  equations 
with  three  unknowns : 


6.90C 


■  32.20C  -  46.00C  =0 

aa  ar  rr 

32.00C  +  66.00C  +  2.00C  =  0.0085 

aa  ar  rr 

32. 0C  -  62.00C  -  1.00C  =  0.0172 

aa  ar  rr 


(4) 


These  equations  give  these  constants: 
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C  =  452.80  x  10  u  (1/psi) 

3d 
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a 

=  6.9  psi 

Act 
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32.0  psi 

Act 

r 

=  -23.0  psi 

Act 
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= 

1.0  psi 

- 

Ae 

=  0.0096 

Ae 

— 

0.0145 

- 

a 

a 

Ae 

=  -0.0044 

Ae 

-0.0027 

r 

r 

Ae 

V 

=  0  (assumed) 

Ae 

V 

= 

0.0085 
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C  =  -94.82  *  10‘6  (1/psi) 

Crr  =  134.18  x  10  °  (1/psi) 

This  procedure  was  used  to  determine  the  anisotropic  elastic  constants 
presented  in  Table  5. 

Comparison  of  Computed  and  Measured 
Pore  Pressure  Parameters 


29.  According  to  the  elastic  theory,  Skempton's  A  parameter  for 
an  anisotropic  material  should  be  related  to  the  elastic  constants  by 
(see  Equation  A15) 


A  = 


C  + 
aa 


2C 


ar 


aa 


4C  + 
ar 


2C 


rr 


(5) 


Thus,  the  theoretical  pore  pressure  response  is  subject  to  verification 
by  a  comparison  of  the  measured  and  predicted  A  parameters.  Unfortu¬ 
nately,  the  procedure  used  to  compute  the  elastic  constants  makes  use 
of  pore  pressure  values  that,  in  turn,  are  used  to  predict  the  A  pa¬ 
rameter.  It  can  be  shown  that  the  A  parameters  computed  from  the  data 
given  in  Table  5  are  numerically  equal  to  the  corresponding  measured  A 
parameters  regardless  of  the  actual  material  behavior.  Therefore,  to 
obtain  an  independent  verification  of  theoretical  concepts,  an  A  param¬ 
eter  must  be  computed  using  data  only  from  drained  loading  increments. 

30.  To  compute  the  elastic  constants,  it  is  necessary  to  load  the 
specimens  in  both  axial  and  radial  directions.  Therefore,  data  from  two 
independent  drained  increments  are  required.  However,  as  discussed  in 
paragraph  26,  the  absolute  stiffness  of  the  soft  shales  was  altered  with 
each  load  cycle  so  that  measurements  from  one  drained  step  were  not  com¬ 
patible  with  those  of  other  drained  steps.  Indeed,  the  value  of  using 
both  drained  and  undrained  increments  to  evaluate  elastic  constants  is 
derived  from  the  ability  to  obtain  three  independent  relationships 
within  a  single  load  step. 

31.  As  an  alternative  to  Equation  5  for  verification  of  the  the¬ 
ory,  the  A  parameter  was  computed  using  the  relative  strains  measured 
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during  the  initial  isotropic  consolidation.  As  shown  in  Appendix  A 
(Equation  A26) , 


*  -  limy 

in  which  m  equals  Ac  /Ae  (measured  during  isotropic  consolidation). 

3  r 

32.  The  advantage  of  Equation  6  is  that  it  does  not  depend  on  the 
absolute  stiffness  of  the  clay  shale.  The  A  values  determined  by 
this  equation  should  be  comparable  to  those  measured  throughout  the  sub¬ 
sequent  undrained  loading  increments  provided  the  ratio  of  stiffness  in 
the  axial  and  radial  directions  is  not  altered. 

33.  A  value  of  m  was  obtained  for  each  shale  by  plotting 
versus  e  (Figure  11).  During  the  initial  consolidation  for  each 

a 

shale,  the  value  m  represents  the  slope.  A  comparison  of  the  calcu¬ 
lated  A  value  with  the  average  measured  value  is  given  in  the  tabula¬ 
tion  below.  Note  that  agreement  was  generally  good  although  the  mea¬ 
sured  A  values  tend  to  be  somewhat  variable  for  the  soft  shales  as 
indicated  by  the  standard  deviation. 


Shale 

Classification 

a 

Ae  “  ra 
r 

Calculated 

A 

A 

Standard 

Deviation 

Bearpaw 

Stiff 

2.5 

0.56 

0.60 

0.05 

Kincaid 

Stiff 

2.5 

0.56 

0.53 

0.07 

Pierre 

Soft 

3.0 

0.60 

0.57 

0.14 

Quivira 

Soft 

4.0 

0.67 

0.70 

0.13 

34.  As  shown  in  Appendix  A  (Equation  A15),  the  maximum  theoreti¬ 
cal  B  value  that  can  be  measured  is  given  by 


C  +  4C  +  2C 

aa _ ar _ rr 

+  C  +  4C  +  2C 
K  aa  ar  rr 

w 


where 


n  =  porosity  of  soil 

K  =  bulk  modulus  of  water 
w 
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RADIAL  STRAIN 


Smaller  values  of  B  might  be  measured  if  the  specimen  is  not  fully 
saturated.  The  theoretical  and  measured  B  values  are  listed  in  the 
tabulation  below.  The  initial  B  values  indicate  that  all  samples  were 
saturated.  Although  the  duration  of  the  tests  was  long  (6  months  to 
1  year) ,  the  specimens  remained  reasonably  saturated  as  indicated  by  the 
final  B  value. 


B 

B 

B 

Theoretical 

Measured 

Measured 

Shale 

n 

Maximum 

Initial 

Final 

Bearpaw 

0.319 

0.991 

1.00 

0.89 

Kincaid 

0.400 

0.997 

1.00 

0.98 

Pierre 

0.385 

0.996 

1.00 

-- 

Quivira 

0.306 

0.996 

1.00 

0.96 

Note:  K 

W 

is  assumed  to  be  314 

,000  psi. 

Failure  Characteristics 

35.  At  the  conclusion  of  the  stress  path  tests,  the  clay  shale 
specimens  were  loaded  to  failure  in  undrained  compression,  then  un¬ 
loaded.  The  failure  stress-strain  curve  in  Figure  12  is  typical  for 
the  stiff  shales.  The  stress-strain  curve  displays  a  slight  break  at  a 
deviator  stress  of  50  psi  but  remains  approximately  linear  up  to  140  psi. 
Figure  13  is  typical  of  the  stress-strain  curve  for  the  soft  shales. 

Jn  contrast  to  the  stiff  shale,  the  soft  shales  did  not  display  a  dis¬ 
tinct  break  in  the  stress-strain  curve  before  failure.  The  failure 
characteristics  of  the  clay  shales  tested  are  summarized  as  follows: 


a  -  a  " 
a  r 


o-o 
a  r 


Shale 

psi 

percent 

psi 

percent 

Initial 

Final 

Bearpaw 

55 

0.15 

160 

1.4 

0.61 

0.0 

Kincaidt 

60 

0.30 

145 

1.5 

0.59 

0.34 

Pierre 

50 

0.60 

86 

4.0 

0.53 

0.32 

Quivira 

70 

0.65 

110 

1.8 

0.70 

0.20 

Note:  Strain  measured  from  beginning  of  final  loading. 
*  At  initiation  of  specimen  dilation. 

**  Based  on  initial  and  failure  pore  pressures, 
t  Not  taken  to  failure  (maximum  values  obtained). 
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UsTRA/N  IMPOSED 
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AXIAL  STRAIN,  e„(lO"4  ) 

Figure  12.  Deviator  stress  versus  axial  strain  for  loading  of 

Bearpaw  shale  to  failure 
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600 

AXIAL  STRAIN,  €a(l0‘4) 

Figure  13.  Deviator  stress  versus  axial  strain  for  loading 
of  Pierre  shale  to  failure 


36.  The  A  parameters  for  both  the  stiff  and  soft  shales  decrease 

in  value  as  failure  approached  as  indicated  in  the  tabulation  above.  A 

stress-path  plot  is  shown  for  the  stiff  shales  in  Figure  14  and  for  the 

soft  shales  in  Figure  15.  Note  that  for  an  isotropic  elastic  material 

the  effective  stress  path  for  undrained  loading  would  be  a  vertical 

straight  line;  thus  the  initial  slope  reflects  the  higher  stiffness 

radially  rather  than  axially  (see  Appendix  A).  Above  a  deviator  stress 

(ct  -  o  )  of  60  psi  the  trend  of  the  curves  reverses  because  of  the 

more  dilatant  nonelastic  behavior  of  the  clay  shales.  Interestingly, 

the  deviator  stress  (a  -  a  )  at  which  the  effective  stress  paths 

a  r 

exhibited  dilatant  behavior  was  nearly  the  same  for  all  four  shales 
even  though  the  ultimate  strengths  were  much  greater  for  the  stiff 
shales.  It  is  important  to  note  that  the  elastic  theory  cannot  be  used 
to  predict  pore  pressure  response  beyond  the  point  where  the  stress  path 
displays  dilative  behavior  (i.e.  has  positive  slope). 


DEVIATOR  STRESS 


PART  V:  DISCUSSION  OF  LABORATORY  RESULTS 


Elastic  Properties 

37.  The  consistency  of  the  computed  elastic  constants  presented 
in  Table  5  suggests  that  they  represent  actual  material  behavior.  The 
greatest  variation  in  properties  appears  to  occur  in  the  calculated 
axial  stiffness  (E^)  of  the  Pierre  shale.  A  detailed  analysis  of  the 
variation  in  stiffness  with  respect  to  load  step  indicates  that  the 
differences  in  stiffness  are  the  result  of  the  inherent  behavior  of  the 
shale  and  not  the  result  of  random  error.  For  example,  the  apparent  uni¬ 
axial  stress-strain  behavior  can  be  constructed  from  the  values  of  E^ 

in  Table  5.  From  the  reconstructed  stress-strain  curve  for  the  Pierre 
shale  (Figure  16),  two  features  are  at  once  apparent.  First,  consider¬ 
able  hardening  occurred  during  step  2  that  made  it  difficult  to  incor¬ 
porate  data  from  the  initial  loadings  into  subsequent  computations  (see 
paragraph  26).  Second,  when  the  specimen  is  subjected  to  a  cycle  of 
loading  and  unloading,  the  shale  exhibits  a  distinct  hysteresis,  which 
results  in  a  different  modulus  for  each  drained  increment.  Note,  how¬ 
ever,  that  for  equivalent  steps  3  and  7  in  Figure  16  the  moduli  values 
are  nearly  equal,  thus  suggesting  that  the  hysteretic  behavior  of  the 
material  is  repeatable  and  indicating  that  the  shales  are  not  perfectly 
elastic.  The  reconstructed  stress-strain  curve  for  the  Bearpaw  shale 
(Figure  17)  shows  that  in  contrast  to  the  softer  Pierre  shale,  the 
Bearpaw  shale  displayed  little  variability  in  the  axial  stiffness  with 
the  exception  of  step  5. 

38.  Because  of  the  loading  program  used,  the  data  obtained  for 
the  Quivira  and  Kincaid  shales  were  not  suitable  for  analysis  of  the 
hysteretic  behavior.  However,  the  plastic  deformation  induced  into 
both  shales  when  subjected  to  a  cycle  of  undrained  loading  indicates  a 
hysteretic  loop  might  have  been  observed  for  drained  increments  of  load¬ 
ing  (similar  to  the  Bearpaw  and  Pierre  shales).  While  it  is  apparent 
that  the  shales  are  not  perfectly  elastic,  the  inelasticity  will  not 
affect  the  pore  pressure  response  provided  the  shale  does  not  become 


Figure  16.  Elastic  behavior  of  Pierre  shale  under  triaxial  loading 
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dilative  (see  paragraph  36).  Further,  the  behavior  of  the  clay  shales 
should  be  amenable  to  analysis  techniques  based  on  the  elastic  analysis 
provided  monotonic  loadings  are  considered.  Significantly,  it  is  evident 
that  much  of  the  variation  in  the  calculated  elastic  properties  is  re¬ 
lated  to  the  actual  material  behavior  and  that  the  testing  methodology 
is  therefore  reliable. 

39.  One  feature  of  the  calculated  elastic  parameters  presented  in 
Table  5  is  the  degree  of  anisotropy  they  exhibit.  A  good  measure  of  an¬ 
isotropy  is  the  constant  r|  ,  which  is  the  ratio  of  radial  stiffness 
to  axial  stiffness  Ej  since  for  an  isotropic  elastic  material  they 
must  be  equal.  While  it  is  impossible  with  the  tests  performed  to  un¬ 
couple  Vj  and  E^/Ej  »  some  indication  of  the  relative  values  of  E^ 
and  E^  can  be  obtained  by  assuming  a  value  of  .  The  value  of  Vj 

reflects  the  coupling  of  strains  in  the  radial  direction  and  is  there¬ 
fore  not  affected  by  the  anisotropy.  It  is,  in  effect,  an  isotropic 
value  (Parry  1979*).  Assuming  =  0.25  ,  the  following  values  of 
E^/E^  are  obtained. 


Shale 

E,/E30  -  »,)« 

Bearpaw 

0.33 

Kincaid 

0.56 

Pierre 

0.53 

Quivira 

0.30 

nt  =  e3/Ej 

2.3 

1.3 

1.4 

2.5 


**  Average  values  obtained  from  Table  5. 
t  Value  of  Vj  assumed  to  be  0.25. 

These  values  suggest  that  the  ratio  of  E^  to  E^  ranges  from  about 
1.3  to  2.5  for  these  four  shales.  It  is  important  to  note  that  if  the 
maximum  possible  value  of  Vj  (0.5)  was  used  in  the  analysis  above,  the 
shales  would  still  appear  highly  anisotropic. 

40.  The  anisotropy  of  the  shales  is  comparable  to  the  value  of 


R.  H.  G.  Parry.  1979.  "Pore  Pressure  in  Clay  Shales,"  Letter  Re¬ 
port,  Geotechnical  Laboratory,  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station,  CE,  Vicksburg,  Miss. 


2.0  determined  by  Parry  (1976)  for  the  Taylor  shale.  Also,  the  ratio  q 
determined  for  these  shales  is  comparable  to  typical  published  values 
for  other  geologic  materials  (Gerrard  1977,  Gerrard  et  al.  1972,  and 
Gibson  1974).  Further,  the  values  of  r|  given  in  the  tabulation  below 
can  be  corroborated  with  the  strain  ratio  data  obtained  from  the 
isotropic  consolidation.  In  Appendix  A  it  is  shown  that 


m(l  -  Vj) 


can  be  calculated: 


'1  - 

1  + 

v2(m  -  2) 

1  assuming 

=0.25  ,  the 

folli 

Shale 

m* 

V2~' 

n 

Bearpaw 

2.5 

0 

1.9 

Kincaid 

2.5 

0.04 

1.8 

Pierre 

3.0 

0.20 

1.9 

Quivira 

4.0 

0.08 

2.6 

(8) 


*  Values  obtained  from  Fig¬ 
ure  1 1 . 

•'~k  Average  values  listed  in 
Table  5. 


The  calculated  values  of  q  appear  to  be  similiar  to  those  obtained 
directly  from  the  tabulation  in  paragraph  39. 


Pore  Pressure  Parameters 


41.  The  most  striking  aspect  of  the  pore  pressure  response  for 
the  four  shales  is  that  in  spite  of  the  large  differences  in  stiffness 
among  the  shales,  they  all  exhibit  A  values  on  the  order  of  0.6  to  0.7. 
These  values  support  the  finding  of  Parry  (1976)  who  suggested  that  the 
magnitude  of  the  pore  pressure  response  exhibited  by  many  clay  shales 
may  be  the  result  of  elastic  anisotropy.  Also,  Parry's  conclusion  is 
strongly  supported  by  the  correspondence  between  the  parameter  m  and 
A  (see  tabulation  in  paragraph  33). 


42.  The  similarity  in  pore  pressure  response  among  such  diverse- 
materials  is  best  explained  by  considering  the  relationship  between 
elastic  properties  and  A  parameters.  The  A  parameter  was  determined 
within  load  ranges  for  which  the  shales  are  essentially  elastic.  From 
Equation  5,  it  is  noted  that  A  is  a  ratio  of  elastic  constants.  In 
terms  of  the  engineering  constants  from  Appendix  A  (Equation  AlOa)  the 
equation  becomes 


1  -  2v„ 


1  -  4v2  +  2 


(1  -  Vj) 


Inasmuch  as  Poisson's  coefficients,  Vl  and  V2  are  within  a  rel¬ 
atively  narrow  range,  A  essentially  depends  on  the  constant  r|  ,  which 
is  the  ratio  of  axial  to  radial  stiffness.  Therefore,  the  many  factors 
that  influence  stiffness  do  not  affect  the  A  parameter.  That  is,  two 
clay  shales  with  drastically  different  stiffness  and  strength  properties 
could  have  the  same  A  parameter. 

43.  The  possible  numerical  value  of  the  A  parameter  is  further 

limited  by  Equation  9.  For  an  isotropic  material,  r|  =  1.0  ,  =  V2  * 

and  A  =  1/3  .  As  the  degree  of  anisotropy  becomes  greater,  so  does  r|  . 
Parry  (1976)  showed  that  as  r|  attains  a  value  of  2.0,  the  A  parameter 
will  be  on  the  order  of  0.55.  For  an  rj  of  5.0,  A  will  be  on  the  or¬ 
der  of  0.80  where  the  maximum  value  of  A  would  be  1.0.  It  would  ap¬ 
pear  that  for  an  intact  (unfractured)  shale  with  significant  anisotropy 

A  would  be  expected  to  fall  within  the  limited  range  of  0.5  to  0.7. 

44.  Note  that  the  above  analysis  assumes  the  specimen  is  oriented 
with  the  highest  stiffness  occurring  in  the  radial  direction.  If  the 
specimen  is  oriented  differently,  the  A  parameter  would  exhibit  differ¬ 
ent  limits.  The  relationship  between  the  orientation  and  the  A  param¬ 
eter  is  given  in  Appendix  A  where  it  is  shown  that  the  specimen  orien¬ 
tation  used  in  this  study  gives  the  maximum  obtainable  A  value. 
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PART  VI:  DETERMINATION  OF  CONSTRUCTION- INDUCED, 

PORE  WATER  PRESSURE 

Analysis  of  the  Embankment  Construction  Problem 

45.  The  pore  water  pressure  induced  in  foundation  materials  can 
be  related  to  total  stress  changes  by  Equation  1.  Therefore,  to  deter¬ 
mine  the  excess  pore  water  pressure  caused  by  embankment  construction, 
it  is  necessary  to  compute  changes  in  the  total  stress  state  associated 
with  excavation  and  fill  operations.  While  a  number  of  analytical  and 
numerical  procedures  are  available  for  stress  analyses,  two  important 
aspects  of  the  problem  must  be  considered.  First,  the  geometric  condi¬ 
tion  should  not  be  overly  simplified.  Embankment  construction  does  not 
proceed  uniformly  across  a  site,  and  piezometers  are  often  located  near 
the  edges  and  corners  of  fill  sections.  An  analysis  based  on  the  as¬ 
sumption  of  plane  stress,  plane  strain,  or  axisymmetric  conditions  does 
not  accurately  account  for  these  edge  and  corner  effects.  Second,  the 
method  should  be  easily  adapted  for  field  use  where  frequent  comparison 
of  observed  and  predicted  piezometer  levels  are  a  part  of  construction 
monitoring.  A  compromise  between  these  two  conflicting  criteria  was 
obtained  by  modelling  the  true  shape  of  the  loaded  area  and  by  simpli¬ 
fying  assumptions  on  the  manner  in  which  the  embankment  transmits  loads 
to  the  foundation. 

46.  The  basic  problem  shown  in  Figure  18a  is  to  determine  the 
change  in  pore  pressure  at  point  A  resulting  from  the  application  of  an 
arbitrary  layer  of  fill.  To  use  the  procedure  developed  for  this  anal¬ 
ysis,  the  problem  was  simplified  by  making  the  following  assumptions 
(Figure  18b): 

a.  Homogeneous  transversely  isotropic  elastic  material 
through  depth  Z  . 

b.  Infinite  loading  surface. 

c.  Nonrigid  fill  layer  acting  as  independent  point  loads. 

d.  Undrained  response  at  point  A. 

e.  Superposition  of  effects  of  all  loads. 
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a.  Basic  problem 


b.  Assumed  problem 


Figure  18.  Simplification  of  embankment  problem 


47.  The  first  assumption  combines  the  existing  fill,  overburden, 
and  rock  foundation  layers  into  an  equivalent  transversely  isotropic  but 
homogeneous  material.  The  second  assumption  gives  the  previously  placed 
fill  layers  the  capability  to  transfer  the  load  as  though  it  extended 
laterally  to  infinity,  while  the  third  assumption  ignores  the  influence 
of  the  rigidity  of  the  newly  placed  fill  layers  on  the  stress  distribu¬ 
tion.  The  fourth  assumption  requires  that  no  pore  water  pressure  dis¬ 
sipation  occur  during  the  load  increment  application  which,  when  com¬ 
bined  with  the  first  assumption,  permits  superpositions  of  effects  of 
all  loads  (fifth  assumption).  While  it  is  possible  to  speculate  on  the 
relative  error  introduced  by  each  assumption,  the  reliability  of  the 
computational  procedure  based  on  these  assumptions  must  ultimately  be 
based  on  comparisons  with  field  observations. 

48.  The  interpretation  of  the  problem  shown  in  Figure  18  reduces 
its  mathematical  formulation  to  two  fundamental  steps.  First,  the  pore 


water  pressure  induced  by  an  individual  point  load  must  be  expressed  as 
a  function  of  its  location  relative  to  point  A.  Second,  the  point  load 
function  must  be  integrated  over  the  loaded  area.  Several  theoretical 
formulas  have  been  developed  to  determine  the  stress  distribution  in  an 
anisotropic  media  for  a  point  load  (Gerrard  and  Wardle  1973).  The 
derivation  of  a  pore  pressure  distribution  function  is  relatively 
straightforward.  The  integration  of  the  distribution  function  is,  how¬ 
ever,  made  difficult  by  complex  fill  shapes  and  in  general  must  be  per¬ 
formed  numerically.  All  pore  pressure  predictions  presented  in  this 
report  were  based  on  the  computer  program  CURLS  developed  to  integrate 
the  distribution  function  over  arbitary  trapezoidal-shaped  fill  areas. 

As  an  aid  to  general  design  computations,  the  program  was  also  used  to 
develop  an  influence  chart  (Figure  19)  for  a  graphical  determination 
of  pore  pressures.  Appendix  B  presents  a  description  of  the  development 
and  use  of  the  influence  chart;  Appendix  C  input  instructions,  a  pro¬ 
gram  listing  for  CURLS,  and  examples  of  its  use. 

Example  Analysis  -  Hillsdale  Dam 

49.  The  construction  and  piezometer  data  of  Hillsdale  Dam  were 
used  to  verify  the  methodology  developed  for  predicting  construction- 
induced  pore  pressures.  The  Hillsdale  Dam  construction  involved  several 
embankment  segments,  built  at  various  rates  over  six  construction 
seasons  (1976-1981).  Further,  the  piezometers  were  located  at  the 
center  line,  edge,  and  corners  of  fill  segments,  providing  data  for  a 
comprehensive  evaluation  of  the  reliability  of  the  pore  pressure  pre¬ 
diction  techniques.  The  embankment  segments  analyzed  are  the  test  berms*, 
cofferdam,  and  the  main  embankment  at  sta  104+00  (Figure  20). 

Description  of  damsite 

50.  The  damsite  is  on  the  Big  Bull  Creek  approximately  35  miles 
southwest  of  Kansas  City  and  5  miles  northwest  of  Paola,  Kansas.  The 
3000-ft  valley  section  of  the  dam  is  centered  between  the  Big  Bull 
Creek  and  Little  Bull  Creek;  their  confluence  lies  immediately  to  the 

south  of  the  dam  (Figure  20a).  The  subdued  upland  relief  necessitated 
*  The  berm  was  part  of  first-phase  construction  and  not  intended  as  a 
test  fill. 
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Figure  19.  Influence  diagram  for  prediction  of  construction-induced 
pore  pressures  in  a  layered  foundation  material 
(A  =  0.7  and  n  =  2.5) 


relatively  long  abutment  sections  extending  7900  ft  to  the  right  of  Big 
Bull  Creek  and  700  ft  to  the  left  of  the  low  bluff  near  Little  Bull 
Creek.  Within  the  valley  section,  the  embankment  obtains  its  maximum 
height  of  75  ft  (Figure  20b,  c,  and  d). 

51.  The  dam  foundation  shown  in  Figure  20b,  c,  and  d  consists 
of  22  to  30  ft  of  alluvial  overburden  underlain  by  sedimentary  lime¬ 
stones,  sandstones,  and  shales  of  the  Pennsylvania-age  Kansas  City 
Group  (Figure  21).  The  alluvium  predominately  consists  of  lean  to  fat 
clays,  although  up  to  3  ft  of  clayey,  gravelly  sand  commonly  covers  the 
bedrock  surface.  Within  the  valley  section,  the  uppermost  rock  units 
consist  of  4  ft  of  Drum  limestone  and  approximately  12  ft  of  Quivira 
shale.  Weathering  in  .<iuoh  of  the  Drum  limestone  has  opened  closely 
spaced,  wavy  shale  partings  and  vertical  joints,  which  in  some  in¬ 
stances  extend  through  the  limestone  to  the  underlying  Quivira  shale. 


•  •••••••• 
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d.  Section  C-C' 

Figure  20.  Plan  view  and  sections  of  Hillsdale  Dam  (modified  from 
U.  S.  Army  Engineer  District,  Kansas  City  1971) 
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SHALE:  soft  to  occasionally  very  soft,  clayey  to  sandy,  platy,  grey  to 
dark  grey,  with  occasional  carbonaceous  partings  and  limestone  nodules, 
upper  half  of  formation  Includes:  sandstone,  moderately  hard,  fine 
grained,  micaceous,  occasionally  calcareous,  thin  bedded,  grey. 


LIMESTONE:  moderately  hard,  dense  to  finely  crystalline,  argillaceous 
and  fo8sillferous ,  thin  to  medium  wavy  bedded,  light  bluish  grey;  2 
shale  units  separate  member  into  3  limestone  units;  shale  is  soft, 
clayey,  platy,  dark  grey. 


SHALE:  soft,  clayey,  platy,  calcareous,  occasional  spherical  phosphatic 
nodules,  grey. 


LIMESTONE:  moderately  hard,  dense,  fosslliferous,  thick  bedded,  light 
grey. 


SHALE,  SANDSTONE,  AND  SILTSTONE:  lateral  as  well  as  vertical  variations 
in  lithology  with  sandstone  and  siltstone  common  in  upper  half  of  forma¬ 
tion;  shale  is  soft  to  moderately  hard,  clayey  to  silty,  platy  to  mas¬ 
sive,  occasionally  calcareous,  dark  grey  to  green  with  pinkish  limestone 
nodules  in  lower  portion  of  formation;  soft  shale  underclay  with  numerous 
sllckensides  occurs  in  middle  of  formation  beneath  thin  coal  seam;  sand¬ 
stone  and  siltstone  is  moderately  hard,  very  fine  grained,  calcareous, 
thin  bedded,  grey;  siltstone  often  interlamlnated  with  shale. 


LIMESTONE:  moderately  hard,  dense  to  very  finely  crystalline,  thin  to 
medium  bedded,  numerous  green  wavy  shale  partings,  fosslliferous,  light 
grey. 


SHALE  AND  SILTSTONE:  shale  is  soft,  clayey  to  silty,  platy,  dark  grey; 
shale  underclay,  with  occasional  very  soft  partings  and  bands,  and  a 
thin  coal  seam  commonly  occur  within  the  upper  4  feet  of  the  member; 
in  lower  portion  of  member  the  shale  is  often  interlamlnated  with 
moderately  hard,  light  grey  siltstone. 


LIMESTONE:  moderately  hard,  thin-bedded,  dense,  argillaceous,  brownish 
grey  limestone  with  green  shaly  partings  and  bands,  varying  to  a  nodular 
limestone  in  a  green  shale  matrix. 


SHALE:  soft  to  moderately  hard,  clayey  to  silty,  platy,  occasional  cal¬ 
careous  and  siltstone  partings,  dary  grey  to  green-grey. 


LIMESTONE:  moderately  hard,  dense  to  very  finely  crystalline,  thin  wavy 
bedding,  light  brownish  grey  with  occasional  light  blue  mottling;  with 
numerous  dark  grey,  soft  shale  partings  to  beds. 


SHALE:  soft  to  moderately  hard,  clayey  to  silty,  platy,  dark  grey  to 
green-grey,  with  occasional  siltstone  partings. 


Bedrock  Un 

t  Thickness 

Parting 

<0.02' 

Band 

0.02  to  0.2' 

Thin  Bed 

0.2  to  0.5' 

Medium  Bed 

0.5  to  1.0' 

Thick  Bed 

1.0  to  2.0' 

Massive 

>2.0' 

Figure  21.  Generalized  stratigraphic  column  for  rock  units  (modified 
from  U.  S.  Army  Engineer  District,  Kansas  City  1971) 
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The  Quivira  shale  consists  of  an  upper  4-ft  section  of  coal  seams  and 
carbonaceous  shale  layers  with  associated  soft  to  very  soft  shale  under¬ 
clays.  The  remainder  of  the  Quivera  shale  consists  of  interlaminated 
soft  shales  and  moderately  hard  siltstones.  Figure  21  describes  the 
rock  units  extending  below  the  Quivira  shale,  which  alternately  consist 
of  limestones  and  shales. 

52.  Treatment  of  the  foundation  prior  to  embankment  construction 
consisted  of  a  grout  curtain,  extended  into  the  rock  units  below  the 
Quivira  shale,  and  a  cutoff  trench,  excavated  through  the  permeable 
overburden  and  Drum  limestone  into  the  Quivera  shale  (Figure  20c  and  d). 
Analyses  of  pore  water  pressure 

53.  The  piezometers  used  in  the  analyses  of  the  berm,  cofferdam, 

and  main  embankment  were  all  located  within  the  Quivira  shale  adjacent 
to  the  main  valley  section  of  the  embankment  between  Little  Bull  Creek 
and  Big  Bull  Creek  (Figure  20a).  The  analyses  were  based  on  the  meth¬ 
odology  described  in  Appendix  B  using  the  computer  program  CURLS  listed 
in  Appendix  C.  From  the  properties  listed  for  Quivira  shale  (see  para¬ 
graphs  33  and  34),  pore  pressure  parameters  A  and  B  were  chosen  to 
be  0.7  and  0.99,  respectively.  The  properties  used  to  compute  the 
stress  distribution  were  based  on  the  assumption  that  the  fill,  over¬ 
burden,  and  rock  collectively  behaved  as  a  transversely  isotropic 
material  with  r)  =  2.5  ,  =  0.2  ,  =  0.1  ,  and  =  0.40Ej  . 

These  values  were  determined  from  published  literature  (Gerrard  1977, 
and  Gibson  1974)  to  be  most  representative  of  layered  systems. 

54.  Test  berm.  The  test  berm  (Figure  20a)  was  constructed 
prior  to  the  main  embankment  during  the  1976  construction  season.  The 
piezometric  levels  were  measured  before  and  during  construction  in  pi¬ 
ezometers  P-94-3  and  P-94-2A  (Figure  22a).  The  construction  sequence 
assumed  for  the  analyses  consisted  of  seven  horizontal  layers  (Fig¬ 
ure  22b).  The  thickness  of  each  layer  was  chosen  to  obtain  the  best 
time  sequence  of  construction  activity  as  inferred  from  fill  inspec¬ 
tion  records. 

55.  Excavation  of  the  cutoff  trench  (Figure  22c)  was  initiated 
during  the  same  construction  season  as  the  berm.  To  evaluate  the 


41 


CUTOFF  TRENCH 


a . 


Plan 


b.  Section  A-A' 


CUTOFF  TRENCH 


fA 


2 


STA  96+00 


P-94-ZA  Q 


Tb- 

SCALE  IN  FEET 
100  0  100  200 


SCALE  IN  FEET 
100  0  100  200 


P-94-3  P-94-2A 


I 


of  test  berm 


influence  of  the  trench,  it  was  assumed  that  the  excavated  material 
could  be  modeled  as  negative  (upward)  loads.  Preliminary  computations 
indicated  that  unloading  during  trench  excavation  and  loading  during 
refilling  of  the  trench  had  virtually  no  influence  on  the  pore  water 
pressures  beneath  any  portion  of  the  berm.  Therefore  the  pore  water 
pressures  induced  by  the  excavation  and  filling  of  the  cutoff  trench 
were  not  included  in  subsequent  analyses  of  the  berm's  construction. 

56.  Relatively  good  agreement  was  obtained  between  predictions 
and  piezometric  measurements  from  P-94-3  at  the  central  portion  of  the 
west  end  of  the  berm  (Figure  23).  However,  the  predicted  values  tended 
to  be  higher  than  those  measured,  particularly  for  measurements  taken 
during  September  1976  when  the  rate  of  loading  was  reduced.  Good  agree¬ 
ment  was  also  obtained  between  predicted  and  observed  trends  for  pi¬ 
ezometer  P-94-2A  located  near  the  edge  of  the  berm  (Figure  24).  How¬ 
ever,  in  contrast  to  the  predictions  for  the  central  portion,  predicted 
piezometric  levels  at  the  berm  edge  were  lower  than  those  measured. 

57.  Cofferdam.  The  cofferdam  was  constructed  in  a  relatively 
short  period  of  time,  1  month  versus  3  months,  and  offered  a  good  ex¬ 
ample  of  undrained  response  during  rapid  loading.  Figure  25  shows  the 
predicted  response  for  piezometer  P-78-1  in  the  center  of  the  cofferdam 
to  be  in  good  agreement  with  the  observation  made  near  the  end  of  con¬ 
struction  in  July  1980.  The  piezometric  level  observed  in  late  Sep¬ 
tember  1980,  nearly  2  months  after  completion  of  the  cofferdam,  was 
about  4  ft  lower  than  the  maximum  predicted  level. 

58.  As  predicted,  the  piezometric  level  adjacent  to  the  coffer¬ 
dam  was  observed  in  piezometer  P-75-2  to  have  virtually  no  immediate  re¬ 
sponse  to  cofferdam  construction  (Figure  26).  However,  based  on  the 
piezometer  reading  in  September  1980,  the  piezometric  level  continued  to 
rise  after  the  end  of  construction. 

59.  Main  embankment.  The  construction  of  the  main  embankment  was 
found  to  be  too  complex  to  develop  a  time-water  level  prediction  (simi¬ 
lar  to  Figure  23)  from  available  field  data.  Therefore,  a  comparison 
was  made  between  the  piezometric  levels  at  the  end  of  construction  and 
the  predicted  distribution  of  piezometric  levels  across  the  axis  of  the 
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Figure  23.  Predicted  and  measured  induced  piezometric  pore  water 
elevations  for  piezometer  P-94-3 
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Figure  24.  Predicted  and  measured  induced  piezometric  pore  water 
elevations  for  piezometer  P-94-2A 
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Figure  25.  Predicted  and  measured  induced  piezometric  pore  water 
elevations  for  piezometer  P-78-1 


CNO  OF 

CONSTRUCTION 

**£DICT£D 

• 

— t  — 

/ 

\ 

M£ASU»£0  ^ 

. 

^ NATURAL  f*OVNDWA  TCK  L£V£l 

JUNE  JULY  AUG  SEPT 

TIME  ,  MONTHS 


Figure  26.  Predicted  and  measured  induced  piezometric  pore  water 
elevations  for  piezometer  P-75-2 
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dam.  Thp  predicted  levels  correspond  to  induced  heads  at  el  845.0  and 
832.0  ft,  which,  respectively,  correspond  to  the  top  and  bottom  of  the 
Quivira  shale.  Table  6  presents  the  specific  predictions  made  for  each 
piezometer  location. 

60.  The  comparison  of  trends  shown  in  Figure  27  reveals  that  near 
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Figure  27.  Predicted  and  measured  induced  piezometric  pore 
water  heads  across  axis  of  dam  (r|  =  2.5) 

the  center  line  of  the  embankment,  the  predicted  piezometric  level  is 
greater  than  that  observed,  while  at  the  edges  of  the  embankment  the 
predicted  is  somewhat  less  than  the  observed.  Therefore,  the  compari¬ 
son  corroborates  the  general  observations  made  for  the  test  berm  and 
cofferdam  that  the  prediction  method  overestimated  the  induced  pore 
water  pressure  near  the  embankment  center  line  but  underestimated  the 
pressure  at  the  embankment  fringe. 
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61.  Comparisons  between  observed  and  predicted  values  of  piezo¬ 
metric  levels  indicate  that  the  proposed  computational  procedure  gives  a 
good  general  prediction  of  construction-induced  pore  water  pressures. 
However,  the  general  correspondence  between  prediction  error  and 
piezometer  location  relative  to  the  embankment  center  line  suggests  the 
error  could  be  attributed  to  some  fundamental  cause.  Four  factors  can 
be  identified  that  could  influence  the  comparison  of  predicted  and  meas¬ 
ured  pore  pressure  changes. 

a.  Fluctuations  in  the  natural  groundwater  level  during 
loading. 

b.  Material  parameters  used  in  analysis. 

c.  Total  stress  computation  procedure. 

d.  Dissipation  of  pore  pressures  during  loading. 

The  potential  error  caused  by  these  factors  and  the  distribution  of  this 
error  with  respect  to  piezometer  location  will  be  discussed  in  the  fol¬ 
lowing  paragraphs. 


Groundwater  Fluctuations 


62.  The  groundwater  fluctuations  observed  prior  to  initiation  of 
construction  were  generally  too  small  to  account  for  the  differences 
between  predicted  and  measured  induced  pore  water  pressures.  However, 
changes  in  the  groundwater  conditions  caused  by  construction  could 
account  for  some  of  the  error.  For  example,  the  excavation  and  refill¬ 
ing  of  the  cutoff  trench  was  observed  to  cause  water  level  fluctuations 
in  piezometers  located  in  the  Quivira  shale.  Based  on  computations, 
these  fluctuations  were  too  great  to  be  attributed  to  the  unloading  and 
loading  associated  with  trench  construction.  Moreover,  concomitant 
fluctuations  were  observed  in  piezometers  located  within  the  Drum  lime¬ 
stone,  which  generally  did  not  respond  to  loading.  It  was  concluded 
that  while  little  seepage  was  observed  to  flow  into  the  trench  excava¬ 
tion,  trench  construction  did  have  sufficient  influence  on  the 
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groundwater  flow  conditions  to  alter  the  piezometric  surface.  Thus, 
the  cutoff  trench  excavation  could  have  been  responsible  for  the  level¬ 
ing  off  of  the  piezometric  levels  observed  in  P-94-3  near  the  end  of 
construction  (Figure  23).  Trends  in  piezometric  levels  shown  in  Fig¬ 
ures  24  through  27  should  not  have  been  similarly  affected. 


Material  Parameters 


63.  The  choice  of  material  parameters  enters  into  the  analysis 
in  two  ways.  First,  the  pore  pressure  parameter  A  must  be  properly 
selected  to  model  not  only  the  intact  material  as  tested  in  the  labo¬ 
ratory  but  also  the  in  situ  conditions  that  include  large-scale  struc¬ 
tural  and  stratigraphic  features.  Second,  the  elastic  properties  used 
in  the  computation  of  total  stress  distribution  must  be  selected  to 
give  the  best  representation  of  the  various  foundation  materials. 
Neither  of  these  factors  can  be  evaluated  solely  from  laboratory  test 
data . 


64.  The  A  parameter  has  an  influence  on  both  the  magnitude  and 
distribution  of  pore  water  pressures.  The  approximate  pore  pressure 
response  induced  by  extensive  fills  is  shown  in  the  following  tabulation 


Au/ytt 


Embankment 

Embankment 

A 

Center 

Edge* 

0.33 

0.90 

0.33 

0.50 

0.92 

0.27 

0.55 

0.93 

0.25 

0.60 

0.94 

0.24 

0.65 

0.95 

0.22 

0.70 

0.96 

0.20 

0.75 

0.97 

0. 19 

1.0 

1.0 

0.11 

Note:  Smallest  lateral  dimension  of  fill 
exceeds  50  times  the  depth  of  the 
piezometer. 

*  Distance  from  fill  is  equal  to  one 
to  two  times  the  depth  of  the 
piezometer. 
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For  materials  with  an  A  parameter  approaching  1.0,  the  inducer}  pore 
water  pressure  Au  is  equal  to  the  applied  vertical  stress.  For  such 
materials,  a  fill  having  a  lateral  dimension  of  three  to  four  times  the 
piezometer  depth  will  produce  pore  water  pressures  equal  to  the  applied 
surface  load  y^t  (Table  Bl) .  For  isotropic  materials  having  A  =  0.33 
the  two  horizontal  stresses  each  have  an  influence  equal  to  that  of  the 
vertical  stress.  For  such  materials,  an  extensive  fill  would  produce  a 
pore  water  pressure  of  only  90  percent  of  the  surface  pressure  ytt  .  i 
serve  that  for  the  range  of  A  values  measured  for  the  Quivira  shale 
(0.55  to  0.7),  the  error  introduced  by  uncertainty  in  the  A  parameter 
is  less  than  5  percent. 

65.  The  error  introduced  in  the  analysis  by  an  incorrect  value 
of  A  is  considerably  greater  for  piezometers  located  near  the  edge  of 
an  embankment  than  it  is  for  an  interior  piezometer.  For  example,  for 
piezometers  located  at  a  distance  from  the  embankment  of  one  to  two 
times  the  piezometer  depth,  the  induced  pore  water  pressure  would  be 
0.3ytt  for  an  A  =  0.33  but  O.ly^t  for  A  =  1  .  Thus,  the  effect 
of  reducing  the  A  parameter  on  the  computed  pore  water  pressure  is  to 
decrease  the  computed  pore  water  pressure  near  the  center  but  increase 
the  computed  pore  water  pressure  near  the  embankment  edge.  Thus,  some 
improvement  in  the  comparison  between  predicted  and  observed  piezo¬ 
metric  levels  could  be  achieved  by  using  a  smaller  A  parameter  in  the 
analysis.  However,  the  magnitude  of  the  improvement  would  be  in¬ 
sufficient  to  explain  the  high  pore  water  pressures  observed  adjacent 
to  the  embankment  (Figure  27).  Moreover,  manipulation  of  the  A  param 
eter  to  improve  the  comparisons  for  embankment  edges  would  create  a  les: 
favorable  comparison  for  the  centrally  located  piezometers  as  indi¬ 
cated  in  Figures  23  and  25. 

66.  The  primary  material  property  that  describes  the  anisotropy 
of  the  foundation  system  is  the  ratio  r)  .  From  the  comparisons  in 
Table  6,  it  is  indicated  that  the  pore  water  pressures  computed  for 
the  isotropic  case,  rj  =  1.0  ,  are  virtually  the  same  as  for  t! 
assumed  anisotropic  case,  r|  =  2.5  .  For  the  extreme  anisotropic 
case,  r|  =  10.0  ,  the  computed  pore  water  pressures  are  about 


20  percent  higher  than  the  isotropic  case  in  the  embankment  centers  and 
several  times  the  isotropic  case  at  the  embankment  edge.  Thus,  unless 
extreme  degrees  of  anisotropy  are  envisioned,  the  induced  pore  water 
pressure  is  not  greatly  affected  by  anisotropy.  Further,  to  improve 
the  comparisons  between  predicted  and  observed  pressures,  a  different 
value  of  r)  would  have  to  be  assumed  for  each  piezometer  location. 

In  view  of  the  correspondence  between  the  location  of  the  piezometer 
within  the  loaded  area  and  the  observed  error,  it  is  unlikely  that  ran¬ 
dom  variations  in  anisotropic  properties  would  fully  account  for  dif¬ 
ferences  between  observed  and  predicted  pressures. 


Computation  Procedure 


67.  The  potential  error  created  by  simplification  of  the  stress 
analysis  is  difficult  to  assess  without  invoking  a  significantly  more 
sophisticated  computational  procedure  to  use  as  a  comparison.  However, 
the  systematic  nature  of  the  error  would  suggest  that  the  simplified 
method  possibly  distorts  the  picture  of  the  true  stress  distribution. 

By  incorporating  the  true  embankment  stiffness  and  the  inhomogeneities 
of  the  layered  foundation,  the  stress  predictions  could  possibly  be 
improved.  However,  since  a  more  sophisticated  analysis  would  also  in¬ 
crease  the  number  of  variables  to  be  determined  (and  the  opportunities 
for  data  manipulation),  it  is  doubtful  that  the  reliability  of  a  more 
sophisticated  analysis  could  be  assessed  with  the  available  data. 

68.  Another  potential  problem  in  the  comparison  of  predicted  and 
observed  piezometric  levels  is  the  assumption  that  the  piezometer  acts 
at  a  distinct  point  within  the  foundation  mass.  In  Figure  28,  the  pore 
water  pressure  varies  by  5  to  10  percent  in  the  vicinity  of  the  piezom¬ 
eter  tip,  an  error  that  exceeds  those  associated  with  the  uncertainties 
in  material  properties.  However,  from  the  predicted  values  shown  in 
Figure  26,  the  error  associated  with  pressure  differences  between  the 
top  and  bottom  of  the  Quivira  shale  would  not  explain  the  observed 
difference  between  computed  and  measured  pore  water  pressures. 
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Figure  28.  Comparison  of  typical  piezometer  installation  with 
computed  induced  head  through  depth  of  Quivira  shale 


Pore  Pressure  Dissipation 


69.  In  view  of  the  pore  pressure  dissipation  that  can  occur  dur¬ 
ing  construction,  the  potential  error  introduced  by  the  assumption  of 
undrained  pore  water  pressure  response  can  be  quite  significant.  The 
magnitude  of  the  error  caused  by  dissipation  cannot  be  fully  assessed 
with  the  available  data.  However,  the  general  influence  of  dissipation 
on  the  pore  pressure  response  was  observed.  For  example,  in  Figures  22 
through  26,  the  greatest  deviation  between  observed  and  predicted  values 
occurred  after  construction  rates  wore  either  reduced  or  the  embankment 
topped  out.  Further,  it  was  generally  observed  for  all  piezometers  that 
piezometric  levels  lower  than  the  predicted  value  tended  to  fall  off  with 
time,  while  levels  higher  than  those  predicted  either  continued  to  rise 


or  remained  constant.  In  many  cases,  the  difference  in  the  two  types 
of  performance  was  quite  dramatic  as  shown  in  Figures  29  and  30.  The 
interpretation  of  these  observations  shown  in  Figure  31  is  that  while 
the  shape  of  the  initial  undrained  pore  pressure  distribution  is  similar 
to  that  of  the  embankment,  with  time  the  consolidation  process  causes  a 
flatter  distribution  to  develop  that  extends  well  beyond  the  limits  of 
the  embankment.  The  interpretation  explains  both  the  distribution  and 
the  magnitude  of  the  observed  error  and  is  consistent  with  all  piezom¬ 
etric  data  obtained  at  the  Hillsdale  damsite. 


C»OJ*0»ATt»  itvu 


a  ic»  ai  oi  ssan  cm  a 


Figure  29.  Piezometric  pore  water  response  after  the  completion 
of  construction  for  piezometer  P-94-3 


Figure  30.  Piezometric  pore  water  response  after  the 
completion  of  construction  for  piezometer  P-104-2 
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re  31.  Conceptional  view  of  pore  pressure  redistribution 
with  dissipation 


PART  VIII:  CONCLUSIONS 


70.  A  combined  laboratory  and  field  investigation  demonstrated 
the  feasibility  of  using  the  theory  of  elasticity  for  undrained  load¬ 
ing  of  transversely  isotropic  materials  to  predict  the  pore  pressures 
induced  in  clay  shale  foundations  by  embankment  construction. 


Laboratory  Investigations 


71.  Based  on  the  four  shales  tested  and  a  previous  investigation 
by  Parry  (1976).  the  following  conclusions  are  made: 

a.  Reliable  data  can  be  obtained  from  long-term  triaxial 
tests  on  compact  clay  shales,  provided  special  pre¬ 
cautions  are  made  to  isolate  air  pressure  from  direct 
contact  with  saturating  and  confining  fluids.  Pore 
pressure  B  values  varied  from  1.0  at  initial  saturation 
to  0.9  or  greater  after  one  year  of  testing. 

b.  Although  shales  are  not  perfectly  elastic,  within  the 
range  of  most  field  loadings,  this  inelasticity  is  of 
secondary  importance,  and  shales  can  be  characterized 
as  elastic  materials  for  prediction  of  construction- 
induced  pore  pressures. 

c.  Clay  shales  are  markedly  anisotropic,  exhibiting  a 
stiffness  1.3  to  2.5  greater  in  the  bedding  plane 
than  normal  planes  (E^/Ej  1.3  to  2.5). 

d.  The  pore  pressure  during  undrained  loading  can  be 
related  to  changes  in  total  stress  by 

Au  =  B(A(J  +  A(Ao  -  Ao  )  ] 
r  a  r 

The  pore  pressure  parameters  A  and  B  can  be  related 
to  the  anisotropic  elastic  constants. 

e.  By  combining  the  pore  pressure  response  for  an  undrained 
loading  segment  with  the  subsequent  drained  stress-strain 
response,  a  computational  procedure  was  developed  to  de¬ 
termine  three  of  the  five  elastic  constants  needed  to 
describe  the  transverse-isotropic  behavior.  The  remain¬ 
ing  constants  can  be  adequately  estimated  from  published 
data . 

f.  Theoretical  values  for  the  pore  pressure  parameters, 


based  on  the  elastic  properties,  compared  favorably  with 
those  measured  in  the  laboratory  tests.  Thus,  the 
theory  of  elasticity  gives  a  reliable  basis  for  estimat 
ing  construction-induced  pore  pressures  in  clay  shale 
foundations . 

g.  The  pore  pressure  parameter  A  essentially  depends  upon 
the  ratio  of  anisotropy,  E^/Ej  >  and  not  magnitude 
of  stiffness.  For  elastic  materials,  A  ranges  between 
1/3  for  isotropic  materials  to  1.0  for  extreme  anisotropy 
Since  the  E^/E^  ratio  of  clay  shales  commonly  falls 

within  a  range  of  1.5  to  3.0,  the  A  parameter  for 
natural  transversely  isotropic  materials  tends  to  fall  in 
the  more  restricted  range  of  0.5  to  0.7. 

h.  Near  failure,  the  shales  display  inelastic  dilatant 
behavior  and  the  A  parameter  approaches  zero.  However, 
because  the  mode  of  failure  in  the  field  is  controlled 

by  the  usually  horizontal  planes  of  isotropy,  the  shear 
failure  induced  in  a  triaxial  specimen  may  not  be  repre¬ 
sentative  of  field  behavior. 


Pore  Pressure  Predictions 


72.  A  simplified  procedure  was  developed  to  compute  the  pore  pres¬ 
sure  induced  in  the  foundation  by  fill  placed  at  the  ground  surface.  A 
comparison  between  predicted  and  observed  piezometric  levels  indicated 
that  the  procedure  gives  good  estimates  of  instantaneous  pore  pressure 
response.  Conclusions  on  the  use  of  the  simplified  procedure  include: 

a.  The  material  parameters  having  the  greatest  influence  on 
the  computed  pore  pressures  are  the  parameter  A  for  the 
shale  and  the  E^/E^  ratio  for  the  combined  foundation 

layers.  It  was  determined  that  for  the  limited  range  of 
numerical  values  typically  observed  for  the  parameters, 
differences  between  observed  and  predicted  pore  pressures 
could  not  be  attributed  to  uncertainties  in  material 
properties . 

b.  Natural  and  construction-related  changes  in  groundwater 
levels  had  a  relatively  small  influence  on  induced 
pore  pressures. 

c.  Redistribution  of  pore  pressures  during  and  after  con¬ 
struction  was  shown  to  have  the  greatest  influence  on 
the  differences  in  observed  and  predicted  pore  pressures. 
Therefore,  to  evaluate  the  pore  pressure  response  of 


clay  shale  foundations  better,  the  analysis  should 
consider  the  effects  of  pore  pressure  dissipation  and 
migration. 

73.  The  inevitable  redistribution  of  pore  pressures  has  important 
practical  implications  in  design  and  construction.  It  is  important  to 
recognize  that  pore  pressure  dissipation  does  not  imply  pore  pressure 
reduction  when  the  two-dimensional  aspects  of  the  problem  are  consid¬ 
ered.  At  and  beyond  the  edges  of  the  embankment,  pressures  may  rise  with 
time,  possibly  to  the  detriment  of  foundation  stability.  Pore  pressure 
measurements  should  be  an  important  part  of  the  construction  monitoring 
program,  especially  in  closure  sections  where  high  pore  pressures  can  be 
developed  before  final  embankment  construction.  Refer  to  EM  1110-2-1908, 
entitled  "Instrumentation  of  Earth  and  Rock-fill  Dams,"  part  1,  para¬ 
graph  4-2,  for  guidance. 
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Incre¬ 

ment 


Aaa  A°r  Au  (10~4)  (1Q~4)  AV/V(10~4) 


28.9  a  =  30.0 
r 


u  =  109.0 


43.0 

-40.0 


a  =68.0 
r 


u  =  112.0 


a  =69.0 
r 


u  =  111.0 


0  =  72.0 

r 


u  =  110.0 


-8.0 

-12.0 


-18.0 


44.8  a  =  73.0 
r 


u  =  109.0 


12.0 

-12.0 


67.4  a  =  68.0 
r 


u  =  109.0 


-13.0 

12.0 


-11.0 


a  =39.0 
r 


u  =  107.0 


-6.0 

-12.0 


-21.0 


Note:  All  pressures  are  in  pounds  per  square  inch. 


•  9 


•  9 


Note:  All  pressures  are  in  pounds  per  square  inch. 
*  No  increment. 


pressures  are  in  pounds  per  square  inch. 


Table  5 

Summary  of  Calculated  Parameters 


Shale 

Incre¬ 

ments* 

c 

aa 

10  ^  •  1/psi 

c 

ar 

10  ^  •  1/psi 

c 

rr 

10  ^  •  1/psi 

Ej,  psi 

v  ** 

2 

E1 

E3 

(1  -  Vj) 

Bearpaw 

3-4 

65.11 

5.06 

14.73 

15,400 

-0.08 

0.23 

5-6 

60.06 

12.46 

25.56 

16,700 

-0.21 

0.43 

7-8 

75.19 

-1.60 

25.56 

13,300 

0.02 

0.34 

9-10 

53.21 

-2.99 

22.45 

18,800 

0.06 

0.42 

11-12 

64.81 

1.08 

16.65 

15,400 

-0.02 

0.26 

13-14 

63.15 

4.59 

18.19 

15,800 

-0.07 

0.29 

Kincaid 

1-20 

99.28 

-9.88 

53.59 

10,100 

0.10 

0.54 

3-20 

82.59 

-1.82 

45.80 

12,100 

0.02 

0.56 

15-20 

75.82 

1.46 

42.64 

13,200 

-0.02 

0.56 

23-20 

72.84 

2.90 

41.25 

13,700 

-0.04 

0.57 

24-20 

81.05 

-1.07 

45.08 

12,300 

0.01 

0.56 

26-20 

89.54 

-5.17 

49.05 

11,200 

0.06 

0.55 

Pierre 

3-4 

452.80 

-94.73 

134. 18 

2,200 

0.21 

0.30 

5-6 

110.61 

-27.71 

92.15 

9,000 

0.25 

0.83 

7-8 

338.48 

-69.01 

144.95 

3,000 

0.20 

0.43 

9-10 

165.78 

-27.01 

80.95 

6,000 

0.16 

0.49 

11-12 

274.77 

-59.54 

107.09 

3,600 

0.22 

0.39 

13-14 

109.71 

-16.62 

77.86 

9,100 

0.15 

0.71 

Quivira 

3-6 

206.40 

-14.97 

42.20 

4,900 

0.07 

0.20 

21-22 

162.40 

-10.11 

67.45 

6,200 

0.06 

0.42 

23-24 

143.13 

-13.78 

38.94 

7,000 

0.10 

0.27 

*  Increments  used  to  calculate  parameters  (see  Appendix  B). 

**  Negative  values  imply  errors  in  computing  nearly  zero  values. 
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Table  6 

Predicted and  Measured  Induced  Piezometric  Pore  Water  Heads 


Location* 

Elevation  of  Fill 
at  Piezometer 

Isotropic 

Predicted 

Anisotropic 

Measured 

Piezometer 

X 

Y  Elevation 

n  =  1.0 

n  =  2.5  n 

=  10.0 

P-104-2** 

2435 

500 

837.0 

— 

2.9 

4.6 

8.3 

14.0 

P-104-3 

2418 

320 

840.0 

905.0 

51.9 

54.1 

67.7 

70.0 

P-104-7 

2425 

80 

841.4 

927.2 

99.5 

100.8 

122.5 

81.0 

P-104-13 

2457 

-350 

836.2 

901.7 

40.3 

42.2 

53.6 

57.0 

P-101-2** 

2135 

-575 

837.2 

— 

1.2 

2.1 

4.4 

7.0 

P-78-1 

-140 

360 

849.1 

902.0 

34.9 

35.0 

— 

31.0 

P-94-3 

1380 

350 

844. 1 

902.0 

46.3 

46.8 

— 

32.0 

P-94-2A** 

1495 

575 

846.0 

— 

3.9 

5.8 

— 

6.5 

P-75-2** 

-460 

610 

851.0 

-- 

0.3 

0.6 

-- 

-- 

Note:  All  values  are  in  feet;  elevations  are  in  feet  above  msl . 

*  Center  of  coordinate  axes  is  located  at  x  =  sta  80+15  and  y  =  dam  center  line. 
**  Piezometers  are  located  at  dam  fringe. 


APPENDIX  A:  PORE  PRESSURE  RESPONSE  OF  TRANSVERSELY 
ISOTROPIC  MATERIALS 


[T 


1.  The  relationship  between  pore  pressure  response  and  the  elas¬ 
tic  compressibility  is  important  in  two  respects.  First,  by  establish¬ 
ing  a  relationship  between  the  elastic  compressibility  and  pore  pressure 
response,  limits  can  be  placed  on  the  pore  pressure  response  parameters 
by  simply  defining  the  expected  range  of  numerical  values  of  the  elastic 
constants.  Second,  to  relate  the  pore  pressure  response  observed  in  a 
triaxial  test  to  the  response  under  more  general  loading  conditions 
expected  in  field  problems,  a  theoretical  correspondence  between  the 
pore  pressure  parameters  and  the  mechanical  properties  of  the  material 
is  needed.  In  addition,  it  is  important  to  establish  the  stress  and 
strain  limits  for  which  the  theory  of  elastic  materials  can  be  used  to 
approximate  pore  pressure  response.  Thus,  it  is  important  to  establish 
relationships  between  pore  pressure  response  under  undrained  conditions 
and  compressibility  characteristics  under  drained  conditions  to  test  the 
theory  experimentally. 


Preliminary  Considerations 

2.  The  stress-strain  response  for  linear  elastic  materials  can  be 
expressed  in  the  general  incremental  form 


{Ae}  =  [C]  {Ao} 


(Al) 


where 

{Ae}  =  strain  increments 
[C]  =  matrix  of  elastic  coefficients 
{Ao}  =  effective  stress  increments 

It  is  important  to  recognize  that  for  transverse  isotropy  the  elements 
in  [C]  depend  on  the  orientation  of  the  reference  axes  relative  to  the 
transverse  plane  of  isotropy.  For  development  of  the  theory  for  pore 


* 


Al 


pressure  response  of  a  transversely  isotropic  material,  it  is  convenient 
to  use  a  reference  coordinate  system  in  which  one  axis  is  orthogonal  to 
the  plane  of  transverse  isotropy.  Thus,  using  the  coordinate  system 
shown  in  Figure  la,  {Ae}  ,  {Aa}  ,  and  [C]  can  be  expressed  as 
follows : 


{A£ }  = 


z 

l/2Ay 


{Ao}  = 


1  /  2Ay 


1/2AV 


[C]  = 


XX 

C 

yx 

c 

yx 

0 

0 

yx 

C 

yy 

C 

yz 

0 

0 

yx 

c 

yz 

C 

yy 

0 

1 

G 

xy 

0 

i 

0 

0 

0 

i 

0 

0 

0 

1 

C 

xy 

0 

0 

0 

0 

,  c 

yz 

,  c 

yy 

,  G 

xy 

,  and 

G  are 

yz 

constants.  The  [C]  matrix  can  also  be  written  in  terms  of  the  engi¬ 
neering  constants  Ej  ,  E^  ,  ,  and  : 


A2 


1 

fl 

V2 

'  E1 

V2 

'  E1 

0 

0 

0 

V2 
'  E1 

1 

E3 

V1 
'  E3 

0 

0 

0 

V2 

'  E1 

V1 
"  E3 

1 

E3 

0 

0 

0 

0 

0 

0 

1 

0 

0 

G13 

0 

0 

0 

0 

2(1  +  Vj) 

0 

E3 

0 

0 

0 

0 

0 

2(1  +  Vj) 

E3  J 

Note  that  when  [C]  is  expressed  in  terms  of  engineering  constants,  it 
can  be  readily  observed  that  five  independent  constants  are  required  to 
describe  the  behavior  of  a  transversely  isotropic  material.  If  [C]  is 
described  with  respect  to  a  coordinate  system  (x',y',z'),  the  coeffi¬ 
cient  terms  will  not  be  the  same  as  those  shown  above.  Rather,  the  co¬ 
efficients  will  be  combinations  of  the  five  engineering  constants  and 
the  cosines  of  transformation  that  relate  (x,y,z)  to  (x',y',z'). 
Importantly,  unless  one  of  the  axes  is  orthogonal  to  the  plane  of 
transverse  isotropy,  the  shear  stresses  ATv 
will  contribute  to  the  normal  strains  Ac 

> 

implying  a  coupling  between  shear  stress  a 

of  the  simplicity  gained  by  eliminating  the  shear-normal  coupling  ef¬ 
fect,  the  (x,y,z)  coordinate  system  has  a  particular  advantage  in  the 
derivation  of  the  pore  pressure  response  equations. 


At 

.  |  ,  At  itf  and 

At  ,  , 

X 

y  x  z 

y  z 

,  Ae  ,  ,  and  Ae  , 

,  thus 

x’ 

’  y  ’  z 

and 

volumetric  strain. 

In  view 

General  Pore  Pressure  Response  Relationships 

3.  The  procedure  used  to  develop  the  pore  pressure  response  rela 
tionship  follows  the  theory  proposed  by  Skempton  (1954)  in  which  the 


'  • 


'  • 


'  • 


'  • 


A3 


volumetric  strain  of  a  material  is  assumed  to  be  equal  to  the  volumetric 
strain  of  the  saturating  fluid  within  the  material's  pore  space.  Skemp- 
ton's  equation  was  derived  for  isotropic  materials  subjected  to  an  incre¬ 
ment  of  total  principal  stress  such  that  Ao^  =  A (i.e.,  as  in  a  tri- 
axial  test).  The  form  of  Skempton's  equation  was  such  that  a  component 
of  pore  pressure  was  caused  by  a  change  in  confining  pressure  A , 
and  a  component  caused  by  the  stress  difference  Aa^  -  A 0^  .  That  is, 


Au  =  B(Acig  +  A(AGj  -  Aa^)) 


where  A  and  B  are  Skempton's  pore  pressure  parameters.  An  equation 
of  form  similar  to  Equation  A2  was  derived  by  Parry  (1976)  for  a  trans¬ 
versely  isotropic  material.  To  illustrate  the  application  of  Skempton's 
approach,  the  derivation  of  the  pore  pressure  equations  for  transversely 
isotropic  material  is  presented  in  the  following  paragraphs. 

4.  From  Equation  Al,  the  elastic  volume  change  is  given  by 


Ae  =  Ae  +  Ae  +  Ae  =  Act  (C  +  2C  )  +  (Act  +  Ao  )(C  +  C  +  C  ) 

v  x  y  z  x  xx  yx  y  z  yx  yz  yy 

where  the  generalized  elastic  constants  C  ,  C  ,  and  C  are  used 

xx  yx  ’  yy 

for  conciseness  of  expression.  Note  that  because  of  the  choice  of  the 

(x,y,z)  coordinate  system,  the  terms  G  ,  G  ,  and  G  do  not 

xy  xz  yz 

appear  in  the  expression  for  Ae 


a  =  (C  +  2C  ) 
xx  yx 

b  =  (C  +  C  +  C  ) 
yx  yz  yy 


to  get 


Ae  =  Act  a  +  (Act  +  Act  )b 
v  x  y  z 


If  Act  =  Act  =  Act  =  Ap  and  At  =  At  =  At  =0  (hydrostatic 
x  y  z  r  xy  yz  xz 

stress  increment)  is  applied,  then  the  volumetric  strain  is  given 
by: 


tc 


Aey  =  Ap(a  +  2b)  =  ApCg 


(A4) 


where 


Cg  =  effective  stress  bulk  compressibility 

=  C  +  4C  +  2(C  +  C  ) 

xx  yx  yy  yz 


Since  Aey  and  p  are  independent  of  how  the  sample  is  oriented,  C£ 

is  an  invariant  quantity.  That  is,  C  +  4C  +  2(C  +  C  )  = 

n  J  xx  yx  yy  yz 

Cx'x'  +  4Cy'x'  +  2(-Cy'y'  +  Cy'z'^  ’  where  (x',y',z')  refer  to  an 
arbitrary  reference  axis.  The  bulk  modulus  of  the  material  can  be 
defined  as 


K  =  !_ 

Ks  C 

s 


(A5) 


5.  For  undrained  loading,  if  the  compression  of  the  individual 
grains  is  ignored,  it  can  be  assumed  that  the  volume  change  of  the  pore 
fluid  is  equal  to  the  total  volume  change  of  the  material.  The  volume 
change  of  the  fluid  can  be  determined  from  the  change  in  pore  pressure 
as 


Ae  = 

v 


Au  n 
K 

w 


(A6) 


where 


Au 

n 

K 

w 


change  in  pore  water  pressure 
porosity  of  soil 
bulk  modulus  of  water 


6.  By  equating  Equations  A3  and  A6  and  invoking  the  definition  of 
effective  stress,  the  following  equation  can  be  written  in  terms  of  total 
stress  and  Au  : 


Ae  = 

v 


Au  n 


w 


=  a(Acrx  -  Au)  +  ^CAcr^  "  Au)  +  ^°x  ~ 


(A7) 


A5 


i 


J 

s 

1 


1 


I 


« 


Then  solving  for  Au 


aAa  +  (Act  +  A a  )b 
a  X  V  2 

Au  = - * - 

+  (a  +  2b) 

w 


(A8) 


|  (a  +  2b)  (Aay  +  Actz)  -t-ajAp^  -  |  (Aay  +  Actz)| 


n 


+  (a  +  2b) 


(A8a) 


w 


Using  the  definition  of  C  (Equation  A4),  Equation  A8a  can  be  written 

s 

as 


Au  = 


+  C 


a_ 

C 


Act 


I 


+  Act 


4 


w 


7.  The  above  equation  is  in  the  form  of  Skempton's  equation 
(Skempton  1954)  for  an  isotropic  elastic  media  in  which  the  pore 
pressure  parameters  A  and  B  are  written  as 


B  = 


5-  +  C 
K  s 
w 


with 


A  =  — 

C 

s 


Au  =  B  (Act  +  Act  )  +  A  [act 
M  y  z  Lx 


I  (*y 


+  toz>]  | 


(A9) 


8.  Note  that  B  is  invariant  and  of  the  same  form  as  for  an 
isotropic  material.  However,  the  A  parameter  is  not  invariant,  and 
Equation  A9  is  valid  only  if  the  coordinate  axes  (x,y,z)  are  specified 
relative  to  the  plane  of  transverse  isotropy  as  indicated  in 
Figure  la. 


A6 


% 


•  V 


9.  The  A  parameter  can  be  written  as 


A  = 


C  +  2C 
xx  yx 


C  +  4C  +  2 (C  +  C  ) 
xx  yx  yy  yz 


(A10) 


By  dividing  the  numerator  and  denominator  by  C  and  defining  the 


ratio 


C  E„ 

xx  3 


C  E, 

yy  i 


=  n 


Equation  A10  can  be  expressed  in  terms  of  the  engineering  constants  as 


1  -  2v„ 


A  = 


1  - 


4v2  + 


2(1  -  Vj) 


(AlOa) 


which  is  identical  to  the  anistropic  A  parameter  derived  by  Parry 
(1976). 

10.  Similarly,  the  B  parameter  can  be  expressed  in  terms  of  the 
engineering  constants  as 


B  = 


1  -  4v2  + 


2(1  -  vx) 


n  E, 


K 


1  -  4v2  + 


2(1  -  vj 


(All) 


which  was  also  derived  by  Parry  (1976).  The  bulk  modulus  can  also  be 
related  to  the  engineering  constants  as 


K  = 


'1 


4v2  + 


2(1  -  v,) 


lT 


(A12) 


11.  For  an  isotropic  material  l)  =  1  ,  V2  =  vi  anc*  ^qua- 
tions  AlOa  and  A12  reduce  to 


A7 


•  •  • 


•  • 


m  9 


which  are  well-known  results  of  elastic  theory. 


Special  Stress  and  Boundary  Conditions 


12.  The  relationships  for  anisotropic  materials  can  be  simplified 
for  many  practical  problems  in  which  special  stress  or  boundary  condi¬ 
tions  exist.  Three  important  special  cases  are  the  triaxial  test, 
drained  isotropic  compression,  and  plane  strain  conditions. 

Triaxial  test 

13.  Triaxial  test  with  (o  =  o  and  a  =  a  =  a  )  For  this 

_  x  a _ y  z  r  . 

test,  only  two  stresses  are  controlled  independently.  Thus,  the  stress 
and  strain  increments  are  related  by  only  three  constants  in  the  fol¬ 
lowing  equation: 


Ae  =  C  Ao  +  2C  Act 
a  aa  a  ar  r 


Ae  =  C  Act  +  C  Act 
r  ar  a  rr  r 


(A13) 


These  constants  are  related  to  the  elastic  properties  as  follows: 


C 

aa 


l 


C 


ar 


(A14) 


c  =  c  +  c  =  — = — - 

rr  xy  yy 

Note  that  only  three  independent  constants  are  involved  in  the  triaxial 

test.  However,  C  ,  C  ,  and  C  are  sufficient  to  determine 
’  aa  ’  ar  ’  rr 

A  and  B  since 


This  results  in 


C  =  C  +  4C  +  2C 
s  aa  ar  rr 


a  =  C  +  2C 
aa  ar 


C  +  2C 
_ aa _ ar 

C  +  4C  +  2C 
aa  ar  rr 


C  +  4C  +  2C 
aa _ ar _ rr 

+  C  +  4C  +  2C 

K  aa  ar  rr 
w 


(A15) 


14.  As  an  aid  to  the  interpretation  of  test  data,  it  is  of  inter¬ 
est  to  relate  the  direction  of  the  effective  stress  path  to  the  A 
parameter.  First,  Equation  A9  must  be  written  in  terms  of  the  effective 
stress  and  pore  pressure  as 


Au  =  B  <o(Aa  +  Act 

)2\  y  z 


Then,  solving  for  Au  ,  the  relationship  obtained  is 


;)  +  A  [“x  -  IK  +  K]*  iU| 


(1  -  B) 


=  IK +  +  *  -  iK +  ^  <ai« 


It  is  immediately  seen  that  if  B  =  1.0  ,  the  effective  stresses  are 
related  by 


Act  +  Act 

\  y  z 


)  *  A[A5x  -  IK  +  ° 


For  the  triaxial  test  with  A(Jx  =  ACTa  and  Act^  =  ACTz  =  ACTf  ,  this  rela¬ 
tionship  reduces  to 


(1  -  A) Act  +  AAct  =  0 
r  a 


A9 


L. 

h 


s 


2 


Ao 


Ao 


A  -  1 


(A17) 


15.  The  stress  path  plots  used  in  this  report  are  based  on  the 

stress  difference  (a  -  a  )  and  mean  stress  (1/3)  (a  +2 a  )  shown  in 

a  r  a  r 

Figure  2.  The  slope  of  the  stress  path  is  given  by 


slope  =  3 


Ao  -  Aa 
_ a _ r 

Aa  +  2Aa 
a  r 


Using  Equation  A17,  the  stress  path  slope  is  related  to  the  A  parameter 
by 

A  -  1 

A 

slope  =  3 


-  1 


A  -  1 
A 

3 

3A  -  1 


+  2 


(A18) 


From  Equation  A18,  it  is  noted  that  for  A  >  1/3  the  slope  is  negative 

and  for  A  =  1/3  (as  for  isotropic  material)  the  slope  is  vertical. 

Further,  since  the  slope  of  the  effective  stress  path  depends  only  on 

A  in  the  triaxial  test,  the  direction  of  the  effective  stress  path  is 

independent  of  the  total  stress  path. 

16.  Triaxial  test  (a  t  a  ,  a  i-  a  and  a  =  o  )  .  If  the 

x  a 


a  f  o 
JL  r 


and  a  =  a  ) 
z  r 


plane  of  transverse  isotropy  is  not  orthogonal  to  the  axis  of  a  triaxial 
specimen,  Equation  A2  is  not  valid.  For  example,  even  though  a 


o^  are  principal  stresses, 


and 


are  not  principal 


and 


and 

r 

strains  since  y  t  0  .  Since  the  specimen  depicted  in  Figure  A1 

3  IT 

will  distort  when  subjected  to  the  principal  stress  increments  Ao 

3 

Aor  ,  the  interpretation  of  the  test  is  greatly  complicated. 

17.  While  distortion  of  the  specimen  limits  the  use  of  a  test  in 
which  the  sample  axis  is  not  orthogonal  to  the  plane  of  isotropy,  such 
a  test  could  still  be  used  to  estimate  the  A  parameter.  Further,  in 
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ORIENTATION 
OF  PLANES 
OF  TRANSVERSE 
ISOTROPY - 


Figure  A1 .  Relationship  between  transverse 
plane  and  specimen  axis 

practice  the  true  orientation  of  the  transverse  isotropic  plane  is  not 
known,  and  it  is  of  interest  to  determine  the  effect  of  its  orientation 
on  the  pore  pressure  response. 

18.  The  influence  of  orientation  on  the  A  parameter  can  be  de¬ 
rived  in  a  manner  similar  to  a  more  general  procedure  proposed  by  Baker 
and  Krizek  (1969).  Considering  the  sample  shown  in  Figure  Al,  Aox  , 

Act  ,  and  Act  can  be  determined  from  the  principal  stresses  Act 
y  z  a 

and  Act^  using  Mohr's  circle  of  transformation 


(A19) 


Act  =  Act 
z  r 


All 


1 


Let  the  A  parameter  at  the  reference  configuration  (given  by  Equa¬ 
tion  A15)  be  denoted  Aq  .  Then  since  B  is  invariant,  Equation  A9 
may  be  written 


B  I  I  [iK  ■"  tor)  '  IK  '  tor)  cos  20  +  tor] 

o  |  IK  *  tor)  +  IK  -  tor)  cos  20 

-  |  |j  (Aaa  ♦  to.)  -  i(toa  -  to.)  cos  20  +  Aor  | 


Rearranging  the  terms 


=  B  |actc  +  ^  [( 1  -  cos  20) 

+  Aq(1  +  3  cos  2e|(Aoa 


which  is  in  the  form  of  Skempton's  Equation  A2.  Hence 


Aq  =  ^£(1  “  cos  20)  +  Aq(1  +  3  cos  20)J 


Note  that  for  an  isotropic  specimen,  Ag  =  Aq  =  1/3  .  The  maximum  and 
minimum  values  of  Ag  can  be  obtained  by  noting  that  they  correspond 
respectively  to  the  maximum  and  minimum  values  of  Au  .  Taking  the 
derivative , 

55  =  J  B(1  -  3A  )(Act  -A a  )  sin  20  (A2 

do  2  oar 


The  maximum  and  minimum  values  are  given  by 


=  o 

90 


which  occurs  when  sin  20  =  0  .  or 


i 


r 


r 

P 


i 


0  =  k  |  ,  k  =  0,  1,  2,. . . 

Hence,  from  Equation  A18,  A  is  a  maximum  at  0=0  and  a  minimum  at 


0  =  90  deg.  The  influence  of 

ing  tabulation: 

0  on 

A  can  be  seen 

A0 

from  the 

follow- 

9>  de8 

A  =  0.3333 

0 

0.50 

0.60 

0.70 

1.0 

15 

0.3333 

0.48 

0.57 

0.66 

0.93 

30 

0.3333 

0.44 

0.50 

0.56 

0.75 

45 

0.3333 

0.38 

0.40 

0.43 

0.50 

60 

0.3333 

0.31 

0.30 

0.29 

0.25 

90 

0.3333 

0.25 

0.20 

0.15 

0 

19.  Provided  the  orientation  of  the  normal  to  the  transverse 
isotropic  plane  is  within  15  deg  to  the  specimen  axis,  the  potential 
error  in  testing  an  improperly  oriented  specimen  is  relatively  small. 
However,  as  0  exceeds  30  deg,  the  potential  error  in  the  determina¬ 
tion  of  A  becomes  quite  significant. 

Drained  hydrostatic 

compression  (Aaa  =  Acr^  =  Aaa  =  Ao^) 

20.  If  a  transversely  isotropic  specimen  is  subjected  to  a  hydro¬ 
static  stress  increment  in  which  Aa  =  Aa  ,  the  ratio  of  m  =  Ae  /Ae 

a  r  a  r 

is  indicative  of  the  degree  of  anisotropy.  For  an  isotropic  specimen 
m  =  1  .  For  material  in  which  E^/E^  >  1  ,  m  >  1  .  The  relationship 
between  m  and  A  is  important  because  it  offers  a  means  of 
obtaining  an  independent  experimental  verification  of  the  relationship 
between  the  material  constants  and  the  pore  pressure  response. 

21.  The  ratio  of  strains  observed  during  drained  isotropic 
compression  can  be  related  to  the  A  parameter  as  follows.  From 
Equation  A13 


•  • 


A13 


Ae  Aa  C  +  2Aa  C 
_ a  _  a  aa _ r  ar 

A£r  Aa  C  +  A a  C 

a  ar  r  rr 


if  Aa  =  Aa 
r  a 


C  +  2C 
aa _ a£ 

C  +  C 
ar  rr 


(A23) 


By  dividing  by  C  and  noting  C  /C 
33  ar  a 


-V£  ,  Equation  A23  can  be 


written  as 


1  -  2v„ 


rr 

C  V2 
aa 


Solving  for  C  /C 

6  rr  aa 


(A23a) 


C  1  -  2v„  +  mV- 
rr  _  2  2 


(A24) 


Similarly,  from  Equation  A15 


1  -  2v„ 


i  -  ^2  +  2  <T 


Substituting  Equation  A24  into  Equation  A25 . 


(A25) 


1  -  2  V 


2 

1  -  2 v„ 


1  -  4v?  +  2 


*  mV2^ 


Therefore , 


m(l  -  2v2) 

m(  1  -  2V2)  +  2(1  -  2v2) 


•  • 


•  9 


A14 


ID  +  2 


(A26) 


Plane  strain  (de  =  0  ,  cr  f  a  ) 
z  y  z 

22.  Plane  strain  conditions  are  often  assumed  to  exist  in  prob¬ 
lems  involving  long  embankments.  In  contrast  to  the  cases  considered 
previously,  the  plane  strain  condition  involves  both  stress  and  strain 
boundary  conditions.  Thus,  to  use  Equation  A8,  A must  be  computed 
from  Act^  and  A using  the  condition  de^  =  0  .  From  Equation  Ala, 
A ct^  can  be  computed  as 


Ao  =  v„r|Ao  +  v.Act 
z  2  x  1  v 


(A27  ) 


Since  Equation  A27  is  written  in  terms  of  effective  stress,  Equation  A16 
must  be  used.  Substituting  Equation  A27  into  Equation  A16  yields 


: B)  =  m 

B  x 


v2n(l  -  A)  +  a]+  \  Act  (1  -  A) (1  +  Vj) 


(A28) 


Define  new  parameters 


Aj  =  \  v2nd  -  A)  +  A 


A2  =  j  0  -  A) (i  +  Vj) 


and  write  Equation  A28  in  terms  of  total  stress 


Au( 1  =  A  Act  +  A  Act  -  (A  +  A  )Au 

B  lx  2y  l  2 


(A29 ) 


Solving  for  Au 


Au  =  B 


A .Ao  +  A„Act 
lx  2  y 


(Aj  +  A2)B  +  (1  -  B) 


(A30) 


Note  that  for  problems  in  which  one  of  the  strains  is  constrained  at 
a  constant  value,  the  pore  pressure  cannot  be  reduced  to  the  simple 


Skempton's  form  (Equation  A9).  However,  if  B  equals  1.0, 
Equation  A30  can  be  reduced  to 


i 


where 


Au  =  Aa  +  A  (A a  -  Act  ) 
y  psv  x  y 


(A31 ) 


Aps  Aj  + 

which  is  in  Skempton's  form  for  B  =  1.0  .  For  isotropic  conditions 
with  A  =  1/3  ,  Vj  =  =  v  <  0.5  ,  and  r)  =  1  ,  note  that 


and 


A,  =  i  (1  t  V) 

A2  =  |  O  +  v) 


As  for  the  A  parameter  (Equation  AlOa),  the  maximum  value  of 
is  1.0. 


APPENDIX  B:  COMPUTATION  OF  PORE  PRESSURES  INDUCED 
BY  EMBANKMENT  LOADS 


1.  The  computation  of  pore  pressures  induced  by  embankment  loads 
is  based  on  the  assumption  that  the  embankment  can  be  replaced  by  an 
equivalent  load  distribution  applied  on  a  level,  homogeneous,  trans¬ 
versely  isotropic  foundation.  Thus,  the  effect  of  the  internal  stiffness 
of  the  embankment  is  ignored  in  the  analysis.  The  equivalent  load 
distribution  is  assumed  to  be  given  by 

P(x,y)  =  Yth(x,y)  dx  dy  (Bl) 


where 

P(x,y)  =  the  load  on  the  infinitesimal  area  dx  dy  at  location 
(x,y)  on  the  loading  surface 

Y  =  total  unit  weight  of  fill  material 

h(x,y)  =  height  of  fill  at  location  (x,y) 

The  stresses  induced  by  a  surface  point  load  can  be  computed  by  the 
appropriate  formula  and  combined  with  Skempton's  equation  to  give  a 
relationship  for  the  pore  pressure  induced  by  a  point  load.  The  pore 
pressure  induced  by  the  embankment  is  thus  given  by 


V 


•  • 


Au 


Yth(x,y)I(x,y) 


dx  dy 


(B2) 


where  I(x,y)  equals  the  influence  factor  derived  from  the  equations 
for  stresses  induced  by  P(x,y)  . 

2.  The  most  commonly  used  equations  for  determining  the  influence 

factor  I(x,y)  are  based  on  the  Boussinesq  theory  (Craig  1978),  which 

gives  for  a  point  at  (x  ,y  ,z  ) 

°  r  o  o  o 


Bl 


2  n 


CT 

z 


3P(x,y) 


2m 


2 


I  <°r  ♦  V  -  ^ 


3r 


z ( 1  -  2v) 


,  2  .  2. 

(r  +  z  ) 


5/2 


(r2  *  *2) 


3/2 


where  from  Figure  B1 

z  =  depth  of  point  (xo,yo,zQ)  below  ground  surface 

r  =  radial  distance  from  P(x,y)  and  (x  ,y  ,z  ) 

o  o  o 


Figure  Bl.  Stresses  in  elastic  half-space  due  to 
point  load  on  surface  (shearing  stress  r  not  shown) 


Writing  Skempton's  equation  in  the  axisymmetric  form 


Au  =  B 


+  Aa  )  +  A  Ao 

r  2 


2  (Aoe 


+  tor>] 


(BA 


the  integrand  for  Equation  B2  can  be  obtained  by  combining  Equations  B1 
B3,  and  B4. 

3.  The  Boussinesq  equation  is  valid  for  isotropic  materials  and 
involves  only  one  material  parameter,  Poisson's  ratio  V  .  Stress 
formulas  for  anisotropic  materials  are  more  complex  and  involve  four 
independent  constants.  A  number  of  formulas  have  been  developed  for 
determining  the  stresses  caused  by  a  point  load  (Gerrard  and  Wardle 
1973).  The  Westergaard  solution  assumes  extreme  anisotropy  in  which 
lateral  strains  (er,£g)  are  not  permitted  (Craig  1978).  Also,  Barden 
(1963)  presented  a  solution  for  orthotropic  materials  in  which  the 
stresses  Aa^  and  A(Jq  were  determined  approximately.  However,  it 
was  found  that  Act^  and  Ao^  had  a  significant  influence  on  the  com¬ 
puted  pore  pressure  and  using  approximate  solutions  introduced  unneces¬ 
sary  uncertainty  into  the  analysis.  Therefore,  the  following  rigorous 
solution  for  the  point  load  problem  presented  by  Lekhnitskii  (1963) 
was  used: 


where 


a  +  c  +  \(a  +  c)  -  4d 
2d 


a  +  c  -  \(a  +  c j  -  4d 
2d 


q]  =  (h  -  as^)(l  -  as*) 
q2  =  (b  -  aS2)(l  -  as*) 


v2(l  +  Vj) 
i  -  nv? 


=  V2n(G13V2  '  V  +  V1G1 

Gn(1  “  ^2) 


*r3 ' (v- *  V2> 

1  -  nv~ 


f  2  2' 

n  -v. 


In  1  -  nv. 


0.4E^  (assumed) 


4.  The  integrand  in  Equation  B2  is  quite  complex  for  fill  shapes 


(described  by  h(x,y))  ,  which  are  of  practical  interest,  and  a  numerical 
integration  scheme  was  used  to  determine  Au  .  To  simplify  the  numerical 
integration  process,  the  function  h(x,y)I(x,y)  was  mapped  into  a  ref¬ 
erence  square  by  using  a  coordinate  transformation  of  the  form 


x  =  £  N  (a,p)X 
i=l  1 


y  =  £  n  (o,p)Y 

i=l  1  1 


where 


X.,Y.  =  coordinates  of  selected  points  on  fill  boundaries 
corresponding  to  s^  ,  t^ 

s^jt^  =  coordinates  of  selected  points  inside  the  reference  square 
a,p  =  coordinates  of  location  within  unit  square 
N^Cotjp)  =  interpolation  factors. 

5.  Using  the  transformation  Equation  B6,  the  integral  B2  can 
be  written  as 


1  1 


Au  =  /  /  BYth'(a,p)I'(or,p)|j|  da  dp 


■1  -1 


where 


=  Jacobian  of  Transformation 


9x  9y  9x  9y 
9s  9t  9t  9s 


Equation  B7  can  be  evaluated  numerically  by  the  quadrature  formula 

k  k 

a-  =  Byt  Z  L  H.H.h' (a. ,p  )I’ (a. ,p.) |j|  . . 

i= 1  i=l  J  J  J  j 


where 


h',I'  =  h(x,y)  and  I(x,y)  evaluated  at  x(a^,Pj)  ,  yCa^p^) 

a^,p^  =  coordinates  of  integration  sampling  points 

=  weighting  functions 

|j|  . .  =  Jacobian  at  a.  ,  p. 

I  I  i.l  i  1 


B5 


Evaluation  of  Equation  B8  requires  only  that  the  fill  height  h  and 
the  influence  factor  I  be  expressed  as  a  function  of  x  and  y  . 

The  influence  factor  can,  of  course,  be  determined  from  combining 
Equation  B4  with  either  Equation  B3  or  B5.  The  function  describing  fill 
height  can  be  expressed  in  the  form  of  an  interpolation  function 
similar  to  the  coordinate  transformation  formulae.  That  is, 


h'(a,0)  =  £  N.(a,0)t 

i=l  1  1 


(B9) 


where 

h'(of,p)  =  fill  height  at  a  ,  p 

t^  =  fill  height  at  in  Equation  B6 

The  sampling  locations  a ^  and  P^  and  the  weighting  function 
and  H.  are  determined  analytically  to  give  the  best  trade-off  between 
accuracy,  economy,  and  simplicity  of  application  (see  for  example, 
Hornbeck  (1975)).  The  general  integration  procedure  described  above  is 
commonly  used  in  the  evaluation  of  stiffness  matrices  in  the  finite 
element  method,  and  a  detailed  description  of  the  procedure  is  given  by 
Zienkiewicz  (1977). 

6.  A  computer  code  was  developed  to  evaluate  Equation  B8  for 
linear  interpolation  factors  (Equations  B6  and  B9),  using  Gauss  inte¬ 
gration  formulae  for  ,  IT  ,  or  ,  and  p^  .  For  a  given  order  of 
integration,  the  Gauss  integration  procedure  specifies  fixed  values  of 
IT  ,  IT  ,  or  ,  and  P^  ,  which  greatly  simplifies  the  computer  code. 
However,  little  advantage  is  gained  from  the  well-known  efficiency  of 
Gauss  integration  formula  since  Equations  B3  and  B5  are  not  well  behaved 
for  the  small  values  of  r/z  typically  found  in  the  embankment  problem. 
The  "shape"  of  the  function  to  be  integrated  is  determined  by  the  loca¬ 
tion  (xo,yo,zo),  and  none  of  the  commonly  used  integration  formulae, 
which  specify  fixed  values  of  or  ,  p^  ,  are  particularly  well  suited 
for  integrating  the  problem  formulated  in  Equation  B7.  Thus,  the  inte¬ 
gration  required  a  large  number  of  sampling  points  or  ,  P^  at  a 
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relatively  large  computational  cost.  Problems  involving  deep  founda¬ 
tion  layers,  relative  to  the  size  of  loaded  areas,  could  be  solved  at 
a  greatly  reduced  cost  since  fewer  sampling  points  would  be  required 
to  evaluate  the  integral  accurately. 

7.  As  an  alternative  to  the  numerical  integration  of  Equation  B2, 
a  graphical  method  based  on  the  influence  chart  technique  was  developed. 
The  influence  chart  consists  of  a  series  of  concentric  circles,  equally 
subdivided  by  radial  lines,  which  are  drawn  to  scale  with  the  center  of 
the  circles  located  over  the  point  for  which  the  pore  pressure  is  to  be 
computed.  The  radii  of  the  concentric  circles  are  determined  so  that 

if  an  infinite  uniform  load  is  placed  over  the  ground  surface,  each 
annular  loaded  region  between  consecutive  circles  makes  an  equal  con¬ 
tribution  to  the  induced  pore  water  pressure.  Thus,  each  influence 
area  created  by  subdividing  the  circles  with  radial  lines  makes  an 
equal  contribution  to  the  induced  pore  pressure.  By  the  principle  of 
superposition,  the  pore  pressure  induced  by  an  irregularly  loaded  area 
can  be  determined  by  summing  the  contribution  of  each  unit  area. 

That  is, 

Au  =  By  I  £  t.  f.  (BIO) 

i=l 

where 

B  =  Skempton's  B  value 
y  =  total  unit  weight  of  fill 
I  =  1/n 

n  =  number  of  subdivided  influence  areas 
t^  =  average  height  of  fill  in  i^*1  influence  area 
f.  =  fraction  of  influence  area  covered  by  fill 

i 

The  influence  chart  for  pore  pressure  is  similar  to  Newmark's  chart  for 
computation  of  vertical  stress  distribution  below  foundations. 

8.  To  illustrate  the  use  of  the  chart,  the  pore  pressure  induced 
by  a  fill  layer  of  the  test  berm  constructed  at  the  Hillsdale  damsite 
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is  computed  for  a  point  34.9  ft  below  the  current  ground  surface.  To 
determine  the  radial  distance  of  each  influence  circle,  a  plot  is  first 
made  of  r/z  versus  Au/y^t  ,  then  r/z  values  corresponding  to 
10  equal  increments  of  Au/y^t  are  determined  from  the  plot  in  Fig¬ 
ure  B2.  The  values  of  Au/ytt  and  r/z  used  to  construct  the  plot 
were  computed  by  the  computer  code  CURLS  described  in  Appendix  C.  The 
radial  distances  used  in  the  scale  drawing  (Figure  B3)  are  computed 
by  multiplying  each  r/z  value  by  34.9  ft.  The  circles  are  further 
subdivided  by  10  equally  spaced  radial  lines  that  give  an  influence 
value  I  of  0.01.  A  circular  broken  line  is  also  drawn  corresponding 
to  Au/y^t  =  0.85  to  facilitate  subdividing  the  larger  influence  areas. 
Using  the  estimated  values  of  f^  and  t^  shown  in  Figure  B3,  the 
pore  pressure  is  computed  as  follows: 

Y,  ft.  =  (70) (1) (17)  +  (1)  (17)  +  (1)  (17)  +  ( 1) ( 17)  +  (1) (17) 
i=l 

+  ( 1 ) ( 1 7 )  +  (1)  (17)  +  (1)  (17)  +  (1) (16)  +  (1) (15) 

+  (1) (16)  +  (1/2) (17)  +  (1/2) (17)  +  (1/2)(17) 

+  (1/2)(17)  +  (1/2)(17)  +  (1/2)07)  +  (1/2)07) 

+  (1/2)03)  +  (1/2)00)  +  (1/2)03)  +  (1/4)(13) 

+  (3/8)(15)  +  (1/4)(12)  +  (1/8) (7)  +  (1/4)(12) 

+  (1/2)07)  +  (1/2)(13)  +  (1/8)  (6)  +  (1/8)  (6) 

+  (1/16)(16)  +  (1/16)04) 

=  1467.6  ft 

For 

B  =  1.0  (assumed) 
yt  =  126.0  (assumed) 

I  =  0.01 

Au  =  (1.0)(0. 01)026. 0)(1467. 6)  =  1850  psf 
Change  in  piezeometric  level  =  =  29.6  ft. 


( 

I! 
R 

i 

E 

R 


» 


r/z 


Figure  B2. 


Au/y^t  values  versus 
problem,  A  =  0.7 


r/z  ratios  used  for  example 
and  r)  =  2.5 


Figure  B3.  Example  problem  showing  the  use  of  the 
influence  chart 
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9.  The  computed  value  of  29.6  ft  compares  well  with  the  observed 

value  of  27.0  ft.  The  value  computed  directly  by  CURLS  was  29.0  ft. 

10.  Data  to  construct  influence  charts  for  other  values  of  A  and 

r|  can  be  computed  using  CURLS  as  described  in  Appendix  C.  However, 

values  of  Au/y^t  for  various  values  of  r/z  and  A  are  tabulated  in 

Table  B1  and  plotted  in  Figure  B4.  Also  shown  in  Table  B1  are  values 

of  Aa  for  various  r/z  ratios.  Thus,  these  data  can  be  used  to 
v 

construct  influence  charts  for  computing  vertical  stresses. 


Table  B1 

Au/y^t  Values  for  Selected  r/z  Ratios  and  A  Parameters 


r/2 

r 

0 

.33 

0, 

.50 

JL 

.55 

_0_ 

.60 

JL 

ml 

\£> 

_0, 

.70 

0. 

.75 

1  . 

,00 

A  a 

V 

0. 

2 

2. 

.  1 

3. 

.  1 

3. 

.3 

3. 

.6 

3. 

,9 

4. 

,2 

4. 

.5 

5. 

9 

5 

.92 

0. 

4 

7. 

.2 

10. 

.3 

11. 

.3 

12, 

.2 

13. 

.  1 

14. 

,  1 

15. 

.0 

19. 

.8 

19 

.77 

0. 

8 

18. 

.6 

26. 

.  1 

28. 

.3 

30. 

.6 

32. 

,8 

35. 

.  1 

37. 

.3 

48. 

.5 

48 

.52 

1. 

6 

35. 

.2 

46. 

.  1 

49. 

.4 

52. 

,7 

56. 

.0 

59. 

.2 

62. 

.5 

78. 

.9 

78 

.94 

3. 

2 

54. 

.0 

64. 

.2 

67. 

.3 

70. 

.3 

73. 

.4 

76. 

.4 

79. 

.5 

94. 

.8 

94 

.80 

6. 

4 

70. 

.5 

77. 

.7 

79. 

.8 

82. 

.0 

84. 

.  1 

86. 

.3 

88. 

.4 

99. 

.  1 

99 

.  12 

12. 

8 

81. 

.  1 

85. 

.8 

87. 

.2 

88. 

.6 

90. 

.0 

91. 

.4 

92. 

.8 

99. 

.9 

99 

.86 

25. 

6 

86. 

.7 

90. 

.0 

91. 

.0 

92, 

.0 

93. 

,0 

94. 

.0 

95. 

.0 

100. 

.0 

99 

.96 

51. 

2 

89. 

.6 

92. 

.2 

93. 

.0 

93. 

.8 

94. 

.5 

95. 

.3 

96. 

.  1 

100. 

.0 

99 

.97 

102. 

4 

91 . 

.  1 

93. 

.3 

94. 

.0 

94, 

.6 

95. 

.3 

96. 

.0 

96. 

.6 

100. 

.0 

99 

.98 

Note:  B  =  1.0  ,  r]  =  2.5  ,  Vj  =  0.2  ,  and  =  0.1  . 


Example: 


Au/y^t  is  in  percent;  O v  is  in  pounds  per  square  foot. 

For  r/z  =  3.2  ,  A  =  0.70 

Au/y^t  =  76.4  percent. 


•  • 


•  • 


•  • 


•  • 


APPENDIX  C:  COMPUTER  PROGRAM  CURLS 


1.  This  appendix  presents  input  instructions  and  computer  code 
for  the  program  CURLS  with  examples  of  its  use.  The  code  is  documented 
with  the  appropriate  equation  numbers  as  found  in  Appendix  B.  The  com¬ 
puter  program  can  be  used  either  for  determination  of  pore-water  pres¬ 
sure  induced  by  an  embankment  load  or  to  produce  the  information 
necessary  for  construction  of  an  influence  chart. 


Input  Instructions 


*  Card  1 _ (8A8) 

Col  5-69:  Title  Any  information  that  is  to  be  printed 

as  a  title  for  the  influence  chart  or 
problem  analysis 

*  Card  2 _ (free  format) 

ITYPE  Type  of  problem: 

0  -  problem  analysis 
1  -  influence  chart 


*  Card  3 _ (free  format) 

XN  Ratio  of  the  horizontal  modulus  to  the 

vertical  modulus  (see  paragraphs  39 
and  54) 

*  Card  4  (free  format) 

A  Value  of  Skempton's  A  parameter 

B  Value  of  Skempton's  B  parameter 

Card  5 _ ( f ree  format ) 

PX  x  coordinate  of  point  where  pore  pres¬ 

sure  prediction  is  to  be  computed  (feet.) 

PY  y  coordinate  of  point  where  pore  pres¬ 

sure  prediction  is  to  be  computed  (feet) 

PZ  Elevation  of  point  where  pore  pressure 

prediction  is  to  be  computed  (feet) 


*  Only  cards  1-4  are  needed  to  obtain  data  for  the  influence  chart. 


Cl 
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Card  6 _ (free  format) 


GAMMA 

Total  unit  weight  of  fill  material 
(pounds  per  cubic  feet) 

NTIMES 

Total  number  of  fill  increments  (integer) 

*  Card  7 

(15) 

NAREA 

Number  of  fill  areas  per  fill  increment 

**  Card 

8 

(free  format) 

XYH (1,1) 

XYH( 1 ,2) 

X  coordinates  of  a  quadrilateral  fill 

XYH(1 ,3) 

XYH( 1 ,4) 

area  input  in  a  clockwise  direction 

**  Card 

9 

(free  format) 

XYH (2,1) 

XYH(2 , 2) 

Y  coordinates  of  a  quadrilateral  fill 

XYH '2, 3) 

XYH (2,4) 

area  input  in  a  clockwise  direction 

**  Card 

10 

(free  format) 

XYH (3,1) 

XYH (3,2) 

Height  of  fill  at  each  corner  of  a 

XYH(3 ,3) 

quadrilateral  fill  area  input  in  a 
clockwise  direction 

XYH(3 ,4) 

**  Card 

11 

(free  format) 

GRADE 

Grade  elevation  before  placement  of  fill 

increment  (feet) 


*  Card  7  is  repeated  for  each  fill  increment  (s.;e  input  data  for 
field  problem). 

**  Cards  8-11  are  repeated  for  each  fill  area  per  fill  increment 
(see  input  data  for  field  problem). 
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Recommended  Material  Parameters 
(Prespecified  in  Program) 

XMU  =  0.2  Poisson's  coefficient  in  the  plane  of 

isotropy  in  a  transversely  isotropic 
material  (line  220) 


XMU1  =  0.1 


Poisson's  coefficient  transverse  to  the 
plane  of  isotropy  in  a  transversely 
isotropic  material  (line  230) 


Computer-Unique  Variables 


The  unit  that  the  computer  reads  input 
data  from  (line  260) 

Sequential  data  line  numbers  (located  in 
each  read  statement) 


Program  Listing 


IOC  ********************************************;!:**#-** 

20C  *  * 

30C  *  CURLS  * 

40C  *  * 

50C  *  PROGRAM  TO  DETERMINE  P.W.P.  INDUCED  IN  * 

60C  *  A  FLAT  LYING*  LAYERED  MATERIAL  DUE  TO  * 

70C  *  EMBANKMENT  LOADING.  PREPARED  BY  ,J,F.  * 

80C  *  PETERS  AT  THE  U.S.A.E.  WATERWAYS  EX-  * 

90  C  *  PER  I  MEN  T  STATION*  VICKSBURG  *  MS.  1981  * 

loot;  *  * 

1  IOC  *******************************************  I:****:** 

120C 

130  DIMENSION  STATE < 50 ) *  WAIT ( 50 ) * XYH ' 3  *  4  )  *  FACT < 4 ) 

140  COMMON  /PROP/GAMMA  *  A  »  B , ISO »  XN*  XMU • XMU1 

150C 

160  CHARACTER  FNA<4>  /*/  • , ■ 33 " > * 33 *  * •  *  V 

170  CHARACTER  TITLE < 10) 

180C 

190C  SAMPLING  POINTS  AND  WEIGHTING  FACTORS  FOR  EQ .  B8 
200  DATA  STATE  /0 . 032380 * 0 . 097005 * 0 . 161222 *0 . 224764 * 0 . 287362 

2108  *  0 . 348756*0. 408686  *  0 . 466903  *  0 . 523161*0. 577225  *  0 . 628867  *  0 . 67787 
2208  0 . 724034  *  0 . 767159  *  0 . 807066  *  0 . 943588  *  0 . 8765720  *  0 . 905879  * 

2308  0,931387*0. 952988  *  0 . 970592 , 0 . 984125 » 0 . 993530  *  0 . 998771  *  26*0 . 0/ 


*  CURLS 

* 

*  PROGRAM  TO  DETERMINE  P.W.P.  INDUCED  IN 

*  A  FLAT  LYING*  LAYERED  MATERIAL  DUE  TO 

*  EMBANKMENT  LOADING.  PREPARED  BY  J.F. 

*  PETERS  AT  THE  U.S.A.E.  WATERWAYS  EX™ 

*  PER I MEN T  STATION*  VICKSBURG *  MS.  1981 

* 

*  *  **  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  %  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  **  *  *  #  t  :!■ ;  *  *  *  1 

DIMENSION  STATE  <  50) *  WAIT (50) *  XYH '3*4) *FACT  <  4 ) 
COMMON  /PROP/GAMMA  *  A  * B  * ISO »  XN  *  XMU * XMU1 

CHARACTER  FNA < 4 )  /*/  ■ , * 33 * > * 33 *  *  *  *  V 

CHARACTER  TITLE < 10) 


250  DATA  WAIT  /0 . 064738 *  0 . 064466  - 0 . 063924  * 0 . 063:1  1 4  *  0 . 062039  * 

260 8  0 . 060704  *  0 . 059 115*0. 057277 , 0 . 055200  *  0 . 052890  *  0 . 050359  *  0 . 0476 1 7 
2708  0.044675, 0.041545, 0.038241  *0.034777  *0.0.?:.  167  *0.027427  *0.023571 
2808  0.019616*0.01 5579  *0.011477*0. 007328  *  0 . 003153  *  26*0 , 0/ 

290C 

300C  PROGRAM  VERSION  FOR  HONEYWELL  SERIES  60  LEVEL  66,  TSS 

3 IOC  RECOMMENDED  MATERIAL  PARAMETERS 

320  XMU=  0.2 

330  XMU1=  0.1 

340C 

35 OC  PROGRAM  READS  FROM  UNIT  12 

360  NF=1 2 

370C 

380C  THE  VARIABLE  LNB  IS  THE  DATA  SET  LINE  NUMBER 
390C 

400C*********** 

410C  FOLLOWING  USED  TO  RUN  FROM  TIME  SHARING 
420C  SET  UP  ACCESS  TO  FILE 

430  PRINT,  "WHAT  IS  YOUR  DATA  FILE  NAME?* 

440  READ  1492*FNA<2) *FNA<3) 

450  1492  FORMAT (2A4) 

460  CALL  ATTACH ( 12*FNA*3*0* ISTAT* ) 

470  CALL  FP ARAM (1,1 32) 

480C*********** 

490C 

500C  FILL  IN  OTHER  HALF  OF  ARRAY 


510  N12=24 

520  N=2#N12 

530  00  2  I = 1 ,  N 1 2 

540  WAIT < I+N12)=WAIT  <  I  > 

550  STATE  < I+N12 )=-STATE  < 1 ) 

560  2  CONTINUE 

570C 

580  9999  FORMAT (//> 

590  1000  FORMAT  < V) 

600C 

610  WRITE <6 ,9999) 

620  READ(NF > 1100)  LNB ,  (TITLE <I)*I=1»8> 

630  WRIT£<6» 1200)  < TITLE< I > , 1=1 »B> 

640  1100  FORMAT ( I4»8A8> 

650  1200  FORMAT  < 1X,66(  ■  ♦  )/,lX,  ,64X ,  ■  *  V  , IX ,  *  *  *  ,  8AB,  -**/,lX, 

6608  ***  ,64X,  **V,1X>66<  *#*  ) > 

670  READ(NF» 1000)  LNB.ITYPE 

680  RE AD (NF, 1000)  LNB  , XN 

690  ISO-  0 

700  XB==  ABS(XN-l) 

710  IFTXB.GT. 0.001)  ISO=  1 

720  REAEi(NF»  1000)  LNB, A, B 

730C 

740C  CHECK  FOR  TYPE  OF  PROBLEM 
750  IF( ITYPE.GT .0)  GO  TO  6 

760C 

770 C  BATA  FOR  FILL 

780  READ ( NF  , 1000 >  LNB,PX,PY*PZ 

790  RE AD < NF » 1 000 )  LNB , GAMMA , NTIMES 

800  WRITE (6,9999) 

810  WRITE (6  * 2000)  PX,PY,PZ 

820  WRITE (6, 3100)  GAMMA, A, B 

830  IF( ISO.EQ.O)  WRITE<6,2100) 

840  IF( ISO.GT.O)  WRITE(6,2200)  XN 

850  WRITE <6, 4000)  NTIMES 

860  WRITE <6, 4001) 

870  GO  TO  7 

880C 

890C  DATA  FOR  INFLUENCE  CHART 

900  6  PX=0.0 

910  PY=0.0 

920  PZ=0.0 

930  NTIMES-  17 

940  GAMMA-1.0 

950  WRITE <6, 9999) 

960  WRITE(6,3000)  A , B 

970  IF ( ISO.EQ.O)  WRITE(6, 2100) 

980  IF < ISO.GT.O)  WRITE ( 6 , 2200)  XN 

990  WRITE <6, 9999) 
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1000  WRITE(6»4050) 

1010  WRITE ( 6  *  4100 ) 

1020C 

1 030C**##***###* 

1040C  BEGIN  LOOP  FOR  TIME 
1050C 

1060  XM=  0.0 

1070  NAREA=1 

1030  7  CONTINUE 

1090C 

1100  TPWP=  0.0 

1110  THEAD=  0.0 

1120  TSV=  0.0 

1130  DO  45  IT=1 f NTIMES 

1140C  CHECK  FOR  TYPE  OF  PROBLEM 

1150  IF ( I  TYPE . GT .  0 )  GO  TO  11 

1160C 

1170C  DATA  FOR  FILL 

1180  READ (NF. 1000)  LNB t NAREA 

1190C 

1200  11  PWP  =0.0 

1210  SV=0.0 

1220  DO  5  K0UNT=1, NAREA 

1230C 

1240C  CHECK  FOR  PROBLEM  TYPE 
1250  IF < ITYPE .EG . 0 )  GO  TO  9 

1260C 

1270C  CALCULATE  DATA  FOR  INFLUENCE  CHART 
1200  IF ( XM .LT.1.0)  XINC=  0.1 


1290 

1300 

1310 

.1320 

1330 

1340 

1350 

1360 

1370 


IF (XM . GE .1.0)  XINC=  XM 
XM=  XM+XINC 

IF<XM.GT. 100.0)  XM*  100.0 
XYH  < 1 » 1 ) =0 . 0 
XYH<1f2)=-0.1736*XM 
XYH  < 1 »  3 ) =0 . 0 
XYH(1»4) =-XYH  <  1  »  2  ) 

XYH ( 2 » 1  ) =0 • 0 
XYH<2,2)=0.9848*XM 


1380 


XYH<2»3)=XM 


1390 

1400C 

14  IOC  1800.0 

1420 

1430 

1440  12 


XYH(2»4)=XYH(2»2) 


GIVES  FACTORS  IN  X  OF  GAMMA*H 
DO  12  1  =  1 » 4 

XYH(3» I ) =1800.0 
CONTINUE 


1450C 


1460C  VARIABLE  XM  EQUALS  R/Z 


1470  GRADE  =1.0 

1480  GO  TO  8 

1490C 


C6 


1500C  READ  DATA  FOR  EACH  AREA  OF  FILL  FOR  TIME  IT 


1510  9  CONTINUE 

1520  DO  10  1=1*3 

1530  READ(NF*1000)  LNB*(XYH(I*J)*J=1*4) 

1540  10  CONTINUE 

1550  READ(NF * 1000)  LNB» GRADE 

1560  8  CONTINUE 

1570C 

1580C  BEGIN  GAUSS  INTEGRATION 
1590C 

1600C  DOUBLE  SUMMATION*  EQ.  B8 
1610  DO  20  K=1 *N 

1620  DO  20  L=1 *N 

1630  XI=STATE (K ) 

1640  ETA=STATE(L) 

1650C  FACT  IS  INTERPOLATION  FACTOR  NI(S.T)!  EQ.  B6 
1660  FACT(1)=<1 .-XI >*(1. FETA) *0.25 

1670  FACT(2)=(1.+XI)*(1+ETA>*0.25 

1680  FACTO  )  =  <1.+XI)*(1.-ETA)*0. 25 

1690  FACT<4)=<1.-XI)*<1.-ETA>*0.25 


X=0. 
Y=0. 
H=0. 
DO  30 


x=x  + 
Y=Y  + 
H=H  + 
CONTINUE 


1  =  1*4 

+  XYH(1*I)#FACT(I> 
+  XYH(2*I)*FACT(I> 
+  XYH(3*I)*FACT(I> 


1670  FACT(2)=(1.+XI)*(1+ETA)*0.25 

1680  FACT(3)=(1.+XI)*(1.-ETA)*0.25 

1690  FACT(4)=(1.-XI)#(1»-ETA)Y0.25 

1700C 

1710C  COMPUTE  X  *  Y  *  AND  H  AT  EACH  INTEGRATION  STATION. 

1720C  USE  EQ.  B6 

1730  X=0. 

1740  Y=0. 

1750  H=0. 

1760  DO  30  1=1*4 

1770  X=X  +  XYH( 1*1) #FACT ( I ) 

1780  Y=Y  +  XYH(2  *  I )*FACT( I ) 

1790  H=H  +  X YH  <  3  *  I ) #FACT ( I ) 

1800  30  CONTINUE 

1810C 

1820C  COMPUTE  P.U.P  DUE  TO  UNIT  LOAD  AT  STATION(K*L> 
1830C 

1840  CALL  PRES (X*Y*H*PX*PY*PZ*  GRADE  *  Dl 

1850C 

1860C  COMPUTE  JACOBIAN  AT  STATION<K*L) 

1870C 

1880  ALLX=XYH< 1 *2)+XYH( 1 *4)-XYH( 1 > 1 )-! 

1890  Al..LY=XYH<2  *  2) FXYH(2*  4  )  ~XYH  ( 2  *  1  >-! 

1900C 

1910  XXI=XYH<  1  *2)+XYH(  1*3)  -XYH<  1*1)  -X' 

1920  YXI=XYH(2»2)+XYH(2*3)-XYH(2* 1 )-X' 

1930  XETA=XYH (1*1) +XYH( 1*2) -XYH (1*3)-) 

1940  YETA=XYH(2*1)FXYH(2*2)-XYH(2*3) -) 

1950C 

1960  XJAKE=XXI*YETA-XETA*YXI 

1970  XJAKE=0.0625*XJAKE 

1980C 


CALL  PRES ( X  * Y*H*PX*PY*  PZ  * GRADE  *  DELU*  SIGZ ) 


C  COMPUTE  JACOBIAN  AT  STATION! K * L > 


ALLX=XYH (1*2) +XYH< 1 *4 )-XYM< 1*1) -XYH (1*3) 
Al..LY=XYH(2 *  2) FXYH(2*  4 ) -XYH ( 2 *  1  )-XYH(2*3) 

XXI=XYH (1*2) IXYH (1*3) -XYH ( 1*1) -XYH (1*4) I  CTA*AI..LX 
YXI=XYH( 2  *  2) +XYH( 2*  3)-XYH( 2  *  1 ) - XYHX  2*4) f ETA* ALLY 
XETA=XYH (1*1) +XYH ( 1*2) -XYH (1*3) -XYH ( 1*4) +XI*ALLX 
YETA=XYH(2*  1  )  +XYH(2*2)-XYH  (2*  3)  -XYH(  2*4)  +XI*ALLY 

XJAKE=XXI*YETA-XETA*YXI 

XJAKE=0.0625*XJAKE 


1990C  Aim  INTO  WEIGHTED  SUM  FOR  INTEGRAL  EVALUATION. 

2000C 

2010  SV=SV  +  SIGZ*XJAKE*WAIT!K>*WAIT(L> 

2020  PWP=PWP  +  DELU*XJAKE*UAIT(K>*WAIT(L) 

2030  20  CONTINUE’ 

2040  5  CONTINUE 

2050  MEAD=0.016*PWP 

2060C  WRITE  OUT  RESULT 

2070  IF ( I TYPE . GT.O)  WRITE! 6. 4200)  XM.PWP.SV 

2080  IF(ITYPE.EQ.O)  WRITE!6.5000)  IT.NAREA.PWP.HEAD.SV 

2090C 

2100C  SUM  VARIABLES  PWP.  HEAD  AND  SV 

2110  TPW P=  TPUP+PWP 

2120  THEAD=  THEAD+HEAD 

2130  TSV=  TSV+SV 

2140  45  CONTINUE 

2150C 

2160C  WRITE  OUT  SUMMATION  OF  PWP.  HEAD  AND  SV 

2170  IF ( ITYPE . EG. 0)  WRITE!6»4300)  TPWP, IHEAD.TCV 

2180C 

2190C#####*###*#*#* 

2200  2000  FORMAT !  IX. “PORE  PRESSURE  INDUCED  AT 

22108  LOCATION  X=  '.F10.2.  *FT.  Y=  '.F10.2.'  FT.'/,5X.'AT  ELEVATION-  ■ 
22208  F8.2» *  FT. ■ ) 

2230  2100  FORMAT! IX. 'MATERIAL  TYPE3  ISOTROPIC) 

2240  2200  FORMAT! IX. 'MATERIAL  TYPE3  ANISOTROPIC.  N=  '.F5.2) 

2250  3000  FORMAT ! IX. *  A  PARAMETER3  1 .F5.3/. IX. *B  PARAMETER3  *.E5.3) 

2260  3100  F0RMAT!1X. 'UNIT  WEIGHT  OF  FILL3  '.F7.1.'  LB/CF’/.IX, 

22708  'A  PARAMETER3  ' . F5 . 3/ . 1 X , • B  PARAMETER3  '.F5,3) 

2280  4000  FORMAT! IX. 'COMPUTATION  FOR*.  15.'  TIME  INCREMENTS'//) 

2290  4001  F0RMAT!3X, 'TIME' »5X. 'NUMBER  OF* »7X. 'PORE' »8X. 'STATIC »4X, 
23008  * VERTICAL '/.IX. 'INCREMENT'. 2X. 'FILL  AREAS  *  »  2X »  * PRESSURE ( PSF ) 

23108  »2X, 'HEAD!FT)'»2X. 'STRESS!  E'SF )*/.lX,9!'-')»2X>10(' -  * ) »  2X . 

23208  13! '-•),2X,8!'-').2X.ll !'-')) 

2330  4050  FORMAT! 12X. 'INFLUENCE  FACTORS ! %  OF  6AMMA*H) '/. 12X . 31 !'-')/ ) 
2340  4100  FORMAT ! 3X. 'R/Z' .5X. 'PORE  PRESSURE* »5X» 'VERTICAL  STRESS'/. 
23508  2X.5! '-' ) .4X.13! ) .5X.15! > ) 

2360  5000  FORMAT !lX»I6»5X»I6»9X»F7.1»6X.F7.2»4X.F7.t) 

2370  4200  FORMAT (1X»F6.1»7X»FA.1» 13X.F7 . 2 ) 

2380  4300  FORMAT !24X»13!*~')»2X»8!'-*)»2X. 11!'-*)/. 16X ."TOTAL3'. 

23908  4X»F8.1,7X»F6«2,3X.F8.1> 

2400  STOP 

2410  END 

2420C 

2430CJF**********  SUBROUTINE  PRES  ***##***##* 

2440C 

2450  SUBROUTINE  PRES!X. Y.H.PX.PY.PZ. GRADE. DELU.SIGZ) 

2460C 

2470C  SUBROUTINE  TO  FIND  P.W.P  INDUCED  IN  A  TRANSVERSELY  ANISOTROPIC 
2480C  CLAY  SHALE  DUE  TO  UNIFORM  LOAD  ON  UNIT  AREA. 


2490C 

2500  COMMON  /PROP/GAMMAr A>B, ISO»XN, XMU,XMU1 

2510  REAL*4  LANDA 

2520C 

2530C 

2540C  COMPUTE  COMMON  FACTORS 
2550C 

2560  SMAL.R=(X- PX)*(X-PX>  +  <Y-PY)*< Y~PY> 

2570  Z=ABS  <  GRADE-PZ ) 

2580  R=SMALR  +  Z*Z 

2590  R=SQRT ( R ) 

2600  SMALR=SQRT ( SMALR ) 

2610  RATIO=SMAL.R/Z 

2620C 

2630C  USE  BOUSSINESQ  EO.  B3 
2640C 

2650  Q=M*GAMMA*0 . 1591549 

2660C 

2670  IF  < ISO. EQ. 1 >  GO  TO  10 

2680  LANDA=3.0*Q*<Z*Z/R**5> 

2690  SIGZ=LANDA*Z 

2700  TAU=LANDA*SMALR 

2710  AVE=SIGZ*  <  3 . 0*SMALR*SMALR-R*R*  <  1 . 0  2 . 0*XMU )  )  /  ( 6 . 0*7*Z  ' 

2720  GO  TO  20 

2730C 

2740C  USE  ANISOTROPIC  EQUATIONS  FROM  LEKHNITSKII 5  EO.  B5 
2750  10  CONTINUE 

2760  XI =- XMU 1  *  ( 1 . 0  f  XMU )  /  <  1 . 0-  XN*XMt.)  1**2) 

2770  X2=< 1 .0-XMU**2)/(XN*(l .0 -XN*XMU1**2> ) 

2780  X3=<0 • 4*XMU1* < 1 . O+XMU) +1.0)/<0.4*(1. 0  -XN*XMU1 **2 > ) 

2790  X0= <  XMU*XN* ( XMU1 -2 . 5  >  +XMU ) / < 1 . 0  XN*XMU1 **2  > 

2800  S3=SQRT(  <  X.1+X3)  **2--4  «0*X2  ) 

2810  S1=SQRT<  <X1+X3+S3)/<2.0*X2>  > 

2820  S2=SQRT ( (X1+X3 -S3)/(2.0*X2> ) 

2830  Ql  =  <  XO -XI *S2**2 ) *  < 1 . 0~X1 *S1**2 ) 

2840  02= ( XO  ~X1*S1**2 )* ( 1 . 0-Xl*S2**2 ) 

2850  X?= ( SMALR*SMALR+S1#S1*Z*Z)**1 . 5 

2860  X6= ( SMALR*SMALR+S2*S2*Z*Z ) **1 .5 

2870C 

2880  SIGZ=(Q*Z/( (S1-S2>*SGRT<X2> ) )*( 1 .0/X9-1 .0/X6) 

2890  AVE=-Q*0.5*Z*< <SQRT(X2)/<Xt*X3-X2> )*( 1 .0/<Rl-82> )*<  S1*S1*Q2/X9 

29008  +S2*S2*Q1/X6 )  +  ( 1 . 0/ ( ( SI -S2)*SQRT  <  X2  > > )*<S1*S1/X9-S2*S2/XA> ) 

2910C 

2920C  LEKHNITSKI I  ASSUMED  COMPRESSION  TO  BE  NEGATIVE i  SO  ........ 

2930  SIGZ=-SIGZ 

2940  AVE=-AVE 

2950  20  CONTINUE 

2960C 

2970C  COMPUTE  PORE  PRESSURE  INCREMENT  FROM  EQ.  B4 
2980  DELU=B*(AVE  +  A* (SIGZ-AVE ) ) 

2990  RETURN 


Example  Problems 

Input  data  for  field  problem 


10  EXAMPLE  PROBLEM  FROM  HILLSBALE!  BAM  -  TEST  BERM 

20  0 
30  2.5 
40  0.7  .999 
50  1380  350  844.1 
60  126.0  7 
70  3 

80  0.0  125.0  185.0  185.0 
90  240.0  490.0  550.0  115.0 
100  3.0  3.0  3.0  3.0 
110  879.0 

120  185.0  185.0  1650.0  1650.0 
130  115.0  550.0  560.0  125.0 
140  3.0  3.0  3.0  3.0 
150  879.0 

160  1650.0  1650.0  1680.0  1760.0 
170  125.0  560.0  500.0  250.0 
180  3.0  3.0  3.0  3.0 
190  879.0 
200  3 

210  10.0  135.0  185.0  185.0 
220  240,0  490.0  535.0  135.0 
230  5.0  5.0  5.0  5.0 
240  882.0 

250  185.0  185.0  1650.0  1650.0 
260  135.0  535.0  540.0  140.0 
270  5.0  5.0  5.0  5.0 
280  882.0 

290  1650.0  1650.0  1680.0  1760.0 
300  140.0  540.0  500.0  250.0 
310  5.0  5.0  5.0  5.0 
320  882.0 
330  3 

340  40.0  160.0  195.0  185.0 
350  240.0  490.0  515.0  150.0 
360  4.0  4.0  4.0  4.0 
370  887.0 

380  185.0  185.0  1650.0  1650.0 
390  150.0  515.0  520.0  160.0 
400  4.0  4.0  4.0  4.0 
410  887.0 

420  1650.0  1650.0  1665.0  1740.0 
430  160.0  520.0  500.0  250.0 
440  4.0  4.0  4.0  4.0 
450  887.0 
460  3 

470  60.0  180.0  185.0  185.0 
480  240.0  490.0  500.0  165.0 
490  5.0  5.0  5.0  5.0 
500  891.0 


510  185.0  185.0  1650. 0  1650.0 
520  165.0  500.0  505.0  175.0 
530  5.0  5.0  5.0  5.0 
540  891.0 

550  1650.0  1650.0  1655.0  1730.0 
560  175.0  505.0  500.0  250.0 
570  5.0  5.0  5.0  5.0 
580  891.0 
590  3 

600  80.0  185.0  185.0  185.0 
610  240.0  445.0  445.0  180.0 
620  2.0  2.0  2.0  2.0 
630  896.0 

640  185.0  185.0  1650.0  1650.0 
650  180.0  445.0  450.0  185.0 
660  2.0  2.0  2.0  2.0 
670  896.0 

680  1650.0  1650.0  1650.0  1715.0 
690  185.0  450.0  450.0  250.0 
700  2.0  2.0  2.0  2.0 
710  896.0 
720  3 

730  90.0  185.0  185.0  185.0 
740  240.0  420.0  420.0  185.0 
750  4.0  4,0  4.0  4.0 
760  898.0 

770  185.0  185.0  1650.0  1650.0 
780  185.0  420.0  430.0  190.0 
790  4.0  4.0  4.0  4.0 
800  898.0 

810  1650.0  1650.0  1650.0  1705.0 
820  190.0  430.0  430.0  250.0 
830  4.0  4,0  4.0  4.0 
840  898.0 
850  6 

860  120.0  185.0  185.0  185.0 
870  240.0  390.0  390.0  240.0 
880  0.0  0.0  0.0  14.0 
890  902.0 

900  185.0  185.0  1650.0  1650.0 
910  240.0  390.0  395.0  250.0 
920  14.0  0.0  0.0  14.0 
930  902.0 

940  1650.0  1650.0  1650.0  1695.0 
950  250  0  395.0  395.0  250.0 
960  14.0  0.0  0.0  0.0 
970  902.0 

980  1650.0  1650.0  1650.0  1695.0 
990  205.0  250.0  250,0  250.0 
1000  0.0  14.0  14.0  0.0 


1010  902.0 

1020  185.0  185.0  1650.0  1650.0 

1030  200.0  240.0  250.0  205.0 

1040  0.0  14.0  14.0  0.0 
1050  902.0 

1060  120.0  185.0  185.0  185.0 

1070  240.0  240.0  240.0  ?00.0 

1080  0.0  14.0  14.0  0.0 
1090  902.0 


* 

Output  data  for  field  problem 


*  ft  IK  #  **  *  *  **  t  #  *  #  #  m  ft  #  #  #  #  #  ft  *  *  *  *  #  *  #  *  A  #  #  *  *  *  #  #  #  *  #  #  #  *  #  #  #  1 ft  ft  ft  ft  ft  *  ft  *  *  *  *  ft  ft  ft  4c  *  * 

ft  ft 

ft  EXAMPLE  PROBLEM  FROM  HILLSDALE  DAM  -  TEST  BERM  * 

*  * 


PORE  PRESSURE  INDUCED  AT  LOCATION  X*  1380. OOFT.  Y=  350.00  FT. 

AT  ELEVATION*  844.10  FT. 

UNIT  WEIGHT  OF  FILL=  126.0  LB/CF 
A  PARAMETER=  0.700 
B  PARAMETER*  0.999 

MATERIAL  TYPE*  ANISOTROPIC*  N=  2.50 
COMPUTATION  FOR  7  TIME  INCREMENTS 


TIME 

NUMBER  OF 

PORE 

STATIC 

VERTICAL 

INCREMENT 

FILL  AREAS 

PRESSURE <PSF) 

HEAD (FT) 

STRESS (PSF 

1 

3 

351.6 

5.63 

398.0 

2 

3 

571.2 

9.14 

651.5 

3 

3 

440.5 

7.05 

509.6 

4 

3 

535.0 

8.56 

627.3 

5 

3 

200.1 

3.20 

241.5 

6 

3 

383.0 

6.13 

468.1 

7 

6 

450.5 

7.21 

541.1 

TOTAL* 


2931.9 


46.91 


3437.1 


Input  data  for  influence  chart 


10  EXAMPLE  PROBLEM  FOR  PRODUCING  INFLUENCE  CHART  FACTORS 

20  1 

30  1.0001 
40  .7  .999 


P 

Output  data  for  Influence  chart 


*####*****###*******#***#*#**##*#***####*#*#**##***************#** 
*  * 

*  EXAMPLE  PROBLEM  FOR  PRODUCING  INFLUENCE  CHART  FACTORS  * 

*  * 
****#****#******##*##*#**#**##****#***#**##*#*# ******************* 


A  PARAMETER*  0.700 
B  PARAMETER*  0.999 
MATERIAL  TYPE*  ISOTROPIC 


l 

I 

l 


INFLUENCE  FACTORS <%  OF  GAMMA*H> 


R/Z 

PORE  PRESSURE 

VERTICAL  STRESS 

0.1 

1.0 

1.47 

0.2 

3.8 

5.68 

0.3 

8.1 

12.07 

0.4 

13.5 

19.87 

0.5 

19.3 

28.33 

0.6 

25.3 

36.81 

0.7 

31.1 

44.87 

0.8 

36.5 

52.23 

0.9 

41.5 

58.78 

1.0 

45.9 

64.49 

2.0 

69.8 

90.98 

4.0 

81,4 

98.54 

8.0 

86.3 

99.78 

16.0 

88  •  6 

99.95 

32.0 

89.8 

99.97 

64.0 

90.3 

99.97 

100.0 

90.5 

99.97 

K 

► 


C13 


•  •  9 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Leavell,  Daniel  A. 

Engineering  properties  of  clay  shales  :  Report  4  : 
laboratory  and  computational  procedures  for  predictions 
of  pore  pressures  in  clay  shale  foundations  /  by 
Daniel  A.  Leavell,  John  F.  Peters,  Frank  C.  Townsend 
(Geotechnical  Laboratory,  U.S.  Army  Engineer  Waterways 
Experiment  Station).  --  Vicksburg,  Miss.  :  The  Station  ; 
Springfield,  Va.  ;  available  from  NTIS,  1982. 

104  p.  in  various  pagings  :  ill.  ;  27  cm.  --  (Technical 
report  ;  S-71-6,  Report  4) 

Cover  title. 

"September  1982." 

"Prepared  for  Office,  Chief  of  Engineers,  U.S.  Army 
under  CWIS  31244." 

Bibliography:  p.  57-58. 

1.  Clay  shales.  2.  Foundations.  3.  Shales. 

4.  Soils--Testing.  I.  Peters,  John  F.  II.  Townsend, 
Frank  C.  III.  United  States.  Army.  Corps  of  Engineers. 
Office  of  the  Chief  of  Engineers.  IV.  U.S.  Army 
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(Card  2) 

Engineer  Waterways  Experiment  Station.  Geotechnical 
Laboratory.  V.  Title  VI.  Series:  Technical  report 
(U.S.  Army  Engineer  Waterways  Experiment  Station)  ; 
S-71-6,  Report  4. 
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